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1.  INTRODUCTION 

The  work  reported  here  is  part  of  a  continuing  program  aimed  at  deter- 
mining the  unsteady  flow  in  a  transonic  compressor  stage.   The  stage  is 
installed  in  the  Turbopropulsion  Laboratories  of  the  Department  of  Aeronautics, 
Naval  Postgraduate  School. 

This  report  has  been  compiled  to  facilitate  use  of  the  data  acquisition 
and  analysis  programs  which  have  been  developed  primarily  for  the  study  of 
unsteady  fluctuating  pressures  on  the  casing  inner  wall.   The  equipment  and 
instrumentation  are  discussed  briefly  in  section  2.   The  data  acquisition 
system  and  programs  are  outlined  in  section  3  and  described  in  detail  in  the 
appendices.   The  post-real  time  data  analysis  programs  are  outlined  with  sample 
results  in  section  4  and  are  described  in  detail  in  the  appendices.   Section 
5  contains  conclusions  and  recommendations  for  further  work. 

2.  EQUIPMENT  AND  INSTRUMENTATION 
2.1  The  Transonic  Compressor 

The  transonic  compressor  test  rig  comprises  an  air  turbine  drive  unit 
and  an  induction  section  which  contains  a  filter,  throttle,  settling  chamber 
and  flow  measuring  nozzle.  The  turbine  drive  unit  supplies  450  HP  at 
30,000  RPM.   The  compressor  is  designed  to  operate  at  30,460  RPM  with  a  rela- 
tive tip  Mach  number  of  1.5.   At  the  design  RPM  and  the  tip  Mach  number,  the 
flow  angle  is  65°  and  the  pressure  ratio  is  1.6  at  a  referred  flow  rate  of 
19  lbm/sec.   The  laboratory  facilities  and  the  test  rig  are  described  in 
detail  by  VAVRA  and  SHREEVE  (1972)  and  VAVRA  (1973). 


2.2  Pressure  Measurement 

Eight  Kulite  CQL-080-25  pressure  transducers  with  natural  frequency  about 
125  kHz  are  mounted  with  their  diaphragms  flush  with  the  inner  case  wall  of 
the  compressor.   Further  details  are  reported  by  PAIGE  (1976).   Table  1  in 
Appendix  A  gives  the  axial  and  circumferential  location  of  the  transducers 
relative  to  transducer  number  K6  which  is  the  furthest  upstream.   The  trans- 
ducers are  used  in  conjunction  with  Datel  Model  201C  instrumentation 
amplifiers  which  have  a  flat  frequency  response  to  100  kHz. 

Each  Kulite  pressure  transducer  is  matched  by  a  pneumatic  static  pressure 
tap  at  the  same  axial  location  in  the  case  wall  (except  in  one  case  -  see  Table 
1  in  Appendix  A)  but  displaced  circumferentially.  Other  pneumatic  static  and 
total  pressure  taps  are  available  upstream.   A  data  recording  system  (VAVRA 
and  SHREEVE,  1972)  is  used  to  record  both  the  steady  pressures  from  the 
pneumatic  taps  and  the  temperature  data.   The  paper  tape  output  from  this 
system  is  processed  using  a  Hewlett  Packard  Model  HP9830A  programmable  calcu- 
lator to  provide  input  data  for  the  measurement  of  fluctuating  pressures  and 
to  establish  the  compressor  operating  point. 

2.3  The  Timing  Disk 

To  enable  synchronization  of  the  sampling  of  the  pressure  transducer 
outputs  with  the  rotation  of  the  rotor,  an  instrumented  timing  disk  is  fitted 
to  the  rotor  shaft.   The  disk  contains  holes  at  intervals  of  one  per  rotor 
blade  and  one  per  rotor  revolution.   Light  sensitive  diodes  and  wave  shaper 
circuits  provide  pulse  trains  to  control  sampling  of  the  pressure  transducers. 
This  system  is  described  in  detail  by  WEST  (1976) . 


3.  DATA  ACQUISITION 
3.1  The  System  Hardware 

Figure  1  is  a  schematic  of  the  data  acquisition  hardware  with  arrows 
indicating  the  flow  of  data  and  control  signals.   The  system  is  under  the 
control  of  the  HP  21MX  computer  which  operates  either  directly  or  through  the 
device  called  "Pacer"  to  control  the  analog-to-digital  (A/D)  converter 
(model  HP5610A)  and  which  transfers  data  to  the  HP9867B  mass  memory  unit  via 
the  HP9830A  calculator. 
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Figure  1.   Schematic  of  the  data  acquisition  system  hardware 


The  peripheral  device  called  the  "pacer"  is  described  in  detail  by 
WEST  (1976)  who  originally  called  it  RPACE.   The  pacer  can  trigger  data 
acquisition  from  a  stationary  transducer  at  any  fixed  point  in  the  rotating 
rotor  frame,  independent  of  the  rotor  speed.   In  effect,  it  divides  the 
circumference  of  the  rotor  into  9  intervals,  each  with  a  circumferential 
length  equal  to  that  of  the  arc  (measured  at  the  wall)  across  two  adjacent 
blade  passages.   Each  of  the  9  intervals  are  subdivided  into  256  equal 
sub-intervals. 

The  pacer  receives  one  per  revolution  and  one  per  blade  input  signals 
from  the  timing  disk  and  performs  two  functions;  it  controls  the  timing  for 
data  acquisition  and  determines  the  speed  of  the  rotor. 

3.2  The  Program  KULITE 

KULITE  is  the  data  acquisition  program  (in  BASIC  language)  for  the 
HP  21MX  computer.   A  flow  diagram,  listing,  variable  assignment  and  notes 
are  presented  in  Appendix  B.   KULITE  can  be  operated  in  three  modes,  viz. 
free- run,  calibration  and  pacer. 

In  free-run  mode  the  A/D  converter  operates  in  mode  4  (see  Reference 
No.  5).   Up  to  1616  samples  of  one  A/D  converter  channel  are  taken  with  a 
frequency  of  10   samples  per  second.   The  sampling  process  is  not  synchronized 
with  the  rotor  rotation. 

In  calibration  mode  the  A/D  converter  operates  in  mode  4  (see  Reference 
No.  5).   It  scans  through  A/D  converter  channel  numbers  1  to  12,  taking 
1616  samples  on  each  channel.   The  average  of  the  1616  samples  is  computed 
before  the  next  channel  is  sampled.   The  scan  is  performed  four  times,  with 
four  different  calibration  pressures  applied  to  the  reference  side  of  each 
Kulite  transducer. 


In  pacer  mode  the  A/D  converter  operates  in  mode  0  (see  Reference  No. 
5).   Sampling  of  a  Kulite  transducer  output  is  synchronized  with  rotor 
rotation  by  means  of  two  pulse  trains  generated  from  light  beams  chopped  by 
the  timing  disk  on  the  compressor  shaft.   One  pulse  train  has  a  frequency  of 
one  per  rotor  revolution  and  the  other  has  a  frequency  of  18  per  rotor 
revolution;  each  pulse  in  the  latter  train  corresponding  to  the  passing  of  a 
blade  past  a  fixed  point.   A  full  description  is  given  by  WEST  (1976). 

In  this  mode  a  pressure  transducer  is  sampled  on  successive  revolutions 
at  a  fixed  point  in  the  rotating  rotor  frame.   Currently,  the  sample  interval 
is  several  revolutions  of  the  rotor.   Changes  in  program  RPACE  would  allow 
samples  to  be  taken  at  intervals  of  one  revolution.   If  the  flow  can  be 
regarded  as  steady  in  the  rotating  rotor  frame  this  technique  enables 
measurement  of  the  wall  pressure  distribution  "carried  around"  by  the  rotor. 
Flow  unsteadiness  in  the  rotating  frame  can  be  averaged  or  the  frequency  con- 
tent of  the  unsteadiness  in  successive  samples  can  be  examined.   In  this 
report  only  averaged  data  from  10  samples  taken  at  each  of  128  points  across 
two  rotor  blade  passages,  is  presented. 

In  all  three  modes  of  operation  the  program  KULITE  transfers  data  from 
the  HP21MX  to  the  HP9830A. 

3.3  The  Program  TRAN4 

TRAN4  is  the  data  acquisition  program  (in  BASIC  language)  for  the  HP9830A 
programmable  calculator.   It  receives  data  from  the  HP21MX  computer,  processes 
it  and  stores  data  on  a  disk  of  the  HP9867B  mass  memory.  A  flow  diagram, 
listing,  variable  assignment  and  notes  are  presented  in  Appendix  C. 


3.4  The  Program  RESET1 

RESET1  initializes  a  record  number  on  the  storage  disk  so  that  at  the 
start  of  a  run  data  can  be  stored  in  file  DATAY1  beginning  at  the  first  record, 
The  program  is  listed  in  Appendix  D. 

4.   DATA  ANALYSIS 
4.1  The  Program 

Off  line  data  analysis  is  at  present  performed  on  the  HP9830A  with  the 
BASIC  language  programs  MAPI,  MAP2,  CONT,  C0NT1,  PLOTSA,  PLOTSB  and  TITIPK. 
These  programs  are  described  in  detail,  with  listings,  flow  diagrams  and 
notes,  in  Appendices  E  through  J. 

MAPI  is  used  to  determine  the  sensitivity  of  the  Kulite  pressure  trans- 
ducers from  data  acquired  with  KULITE  in  the  calibration  mode.  In  addition, 
MAPI  is  used  to  convert  the  voltages  sampled  at  the  pressure  transducer  out- 
puts to  pressure  coefficients. 

MAP2  is  used  to  convert  the  8  x  128  array  of  measured  pressure  coeffi- 
cients to  a  29  x  128  array  through  quadratic  interpolation  in  the  axial 
direction.   The  program  was  written  to  reduce  the  effects  of  the  course  trans- 
ducer spacing  in  the  axial  direction.   However,  care  must  be  exercised  when 
it  is  used  to  interpolated  across  discontinuities  such  as  shock  waves  and 
rotor  blades.   Linear  interpolation  is  available  through  use  of  the  program 
CONT  or  C0NT1  to  plot  contours  of  casing  wall  pressures. 

CONT  is  used  to  plot  contours  of  constant  casing  wall  pressure  (in  the 
frame  of  the  rotor)  from  an  array  of  pressure  coefficients.   (i.e.  it  pro- 
duces a  wall  pressure  "map") .  The  program  will  accept  any  general  rectangular 
array  provided  that  the  spacing  in  each  direction  is  uniform.  This  latter 


requirement  restricts  it's  use  in  this  application  to  arrays  obtained  from  MAP2. 

C0NT1  is  used  to  plot  contours  of  constant  casing  wall  pressure  when  the 
array  of  measured  pressure  coefficients  contains  nonuniform  spacing  in  the 
axial  direction.   Nonuniform  spacing  results  in  this  application  from  the 
axial  location  of  the  Kulite  pressure  transducers. 

PLOTSA  is  used  to  plot  (on  the  HP9862A  plotter)  the  uncalibrated  pressure 
distribution  (in  volts)  across  a  blade  pair  for  a  given  Kulite  transducer. 
The  input  data  is  that  originating  from  pacer  mode  of  operation.   The  program 
is  also  used  to  plot  the  output  of  the  one  per  blade  signal  from  the  timing 
disk. 

PLOTSB  is  used  to  plot  (on  the  HP9862A  plotter)  the  uncalibrated  free- 
run  data  (in  volts)  from  a  given  transducer  against  circumferential  distance. 

TITIPK  is  used  to  superimpose  the  blade  tip  profiles  on  the  wall 
pressure  maps. 

4.2  Sample  Results 

The  results  presented  here  are  intended  only  to  illustrate  the  capabili- 
ties of  the  programs.   Comprehensive  results  will  be  given  in  a  subsequent 
report. 

Figure  2  is  a  plot  versus  circumferential  distance  of  the  average 
pressure  in  the  frame  of  the  rotor  across  an  arc  of  the  casing  wall  equivalent 
to  two  blade  passages.   The  pressure  coefficient  is  defined  in  Appendix  E. 
The  plot  was  made  with  PLOTSA  using  data  acquired  in  the  pacer  mode  of  operatior 
The  precise  location  of  each  distribution  relative  to  the  rotor  blades  is  not 
defined  here.   The  locations  are  known  approximately  from  the  blade  pair 
number  specified  in  the  acquisition  program.   They  are  located  precisely  from 
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Figure  2.   Waveshapes  of  unsteady   pressure  distributions  across  two  blade 
passages  (uncalibrated  and  with  arbitrary  offsets)  for  the  KuUte 
transducers.   Data  taken  in  pacer  mode.   50%  design  speed;  throttled 
to  near  surge.   8.7  Ibm/sec  referred  flow  rate.   Pressure  ratio  = 
1.155:1.   Blade  pair  #2. 
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Figure  3.   Typical  record  of  raw  one  per  blade  signal  taken  in  pacer 
mode.   Blade  passage  frequency  is  4.55  kHz. 
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Figure  4. 


Typical  record  of  data  taken  in  free-run  mode  from  Kulite  number  K10. 
The  record  is  comprised  of  1616  samples  taken  at  a  frequency  of 
100  kHz.   Each  cycle  is  due  to  a  blade  passage  with  a  blade  passing 
frequency  of  4.55  kHz.  Compressor  operating  conditions  as  in  figure  2 
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Figure  5.   Contours  of  constant  pressure  coefficient  ACp  plotted 
by  C0NT1.   (See  Appendix  E  for  definition  of  AC  ) . 
2.   Compressor  operating  conditions  as  in  figuri  2. 
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Figure  6.   Smoothed  contours  of  constant  pressure  coefficient  obtained  from  Figi 
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Figure  7 


Contours  of  constant  pressure  coefficient  plotted  by  CONT  with 
25x255  array.   Blade  pair  number  2.   Compressor  operating 
conditions  as  in  Figure  2. 
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Figure  8.   Contours  of  constant  pressure  coefficient  plotted  by 

CONT  with  7x64  array  as  subset  of  array  used  for  Figure  7 
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Figure  9.   Contours  of  constant  pressure  coefficient  plotted  by  CONT  with 

8x128  array.   Note  there  is  a  small  error  in  blade  location.   Blade 
pair  number  2.   Sixty  percent  design  speed.   Throttled  to  near  surge 
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a  knowledge  of  the  orientation  of  the  timing  disk  relative  to  the  rotor  in 
conjunction  with  the  phase  relationship  between  the  one  per  blade  timing 
signal  and  the  pressure  distributions  from  the  transducers.   A  typical  record 
of  the  one  per  blade  signal,  taken  directly  from  the  photo-diode  output  in 
the  pacer  mode  is  plotted  in  figure  3. 

The  distributions  shown  are  of  output  voltage  from  the  Kulite  transducers 
which  at  this  stage  have  not  been  scaled  to  take  account  of  the  tranducer 
calibrations.   The  rapid  changes  in  the  signals,  in  particular  those  from 
transducers  K9,  K10  and  Kll,  are  due  to  the  passage  of  a  rotor  blade  across  a 
transducer  and  provide  one  means  of  estimating  blade  location.   The  distribu- 
tions have  been  plotted  by  linear  interpolation  between  the  128  points  across 
the  blade  pair.   Each  of  the  128  points  has  been  obtained  from  the  average  of 
ten  samples,  one  sample  being  taken  approximately  each  tenth  revolution  of 
the  rotor. 

Figure  4  is  a  typical  record  of  data  taken  in  the  free-run  mode  from  a 
particular  transducer.   The  fundamental  frequency  is  that  of  the  blade 
passage  past  the  transducer  (4.55kHz).   The  plot  is  a  linear  interpolation 
between  1616  samples  taken  at  a  frequency  of  100  kHz. 

Figure  5  is  a  representative  map  of  the  contours  of  constant  pressure 
coefficient  with  respect  to  an  upstream  reference  static  pneumatic  pressure. 
It  was  generated  with  the  program  C0NT1  from  the  8  x  128  array  of  measurements 
which  were  taken  in  the  pacer  mode  and  calibrated  with  the  program  MAPI.   The 
map  across  two  adjacent  rotor  blade  passages  is  thus  that  of  the  mean  wall 
pressures  "carried  around"  by  the  rotor.   The  method  of  calculating  the 
the  pressure  coefficients  from  the  Kulite  transducer  outputs  and  the  measure- 
ments from  the  pneumatic  taps  is  given  in  Appendix  E. 
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The  location  of  the  blade  tip  in  figure  5  was  determined  from  the  circum- 
ferential location  of  the  rapid  change  in  pressure  distribution  associated 
with  transducer  K9  (figure  2).   The  transducer  K9  has  an  axial  location 
such  that  it  is  crossed  by  the  leading  edge  of  rotor  blades.   A  more  precise 
location  of  blades  can  be  obtained  from  measurement  of  the  circumferential 
location  of  the  one  per  blade  raw  signal  (figure  3)  relative  to  the  pressure 
distributions  (figure  2)  and  from  knowledge  of  the  location  of  the  one  per 
blade  holes  in  the  timing  disk  relative  to  the  rotor. 

Smoothing  of  the  Wall  Pressure  Maps 

It  is  clear  in  figure  5  that  linear  interpolation  between  pressure 
coefficients  in  the  axial  direction  across  a  blade  gives  incorrect  contours 
along  the  blade  between  the  transducer  locations.   This  difficulty  can  only 
be  overcome  satisfactorily  by  using  more  Kulite  transducers  to  provide  a 
finer  mesh  than  the  rather  coarse  one  provided  by  eight  transducers.   In  this 
study  the  contours  in  the  vicinity  of  the  blades  have  been  smoothed  graphically 
by  connecting  with  smooth  curves  those  points  in  the  map  at  which  the  pressure 
coefficient  is  the  same  and  which  lie  on  the  lines  scanned  by  the  transducers. 
At  operating  conditions  which  give  rise  to  shock  waves  it  is  possible  that  a 
similar  smoothing  procedure  will  be  required.   In  the  long  run  there  is  a  need 
for  a  numerical  interpolation  technique  which  avoids  interpolation  across 
blades. 

Figure  7  is  a  wall  pressure  map  produced  by  CONT  for  similar  operating 
conditions  to  those  used  to  obtain  figure  5.   However,  pressures  were 
measured  at  255  points  in  the  circumferential  direction  and  MAP2  was  used  to 
generate  pressure  coefficients  at  25  equally  spaced  stations  in  the  axial 
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direction  between  transducers  K6  and  K12.   It  is  clear  that  quadratic  inter- 
polation does  not  avoid  the  problems  of  contour  distortion  by  interpolation 
across  a  blade.   In  view  of  the  very  slow  execution  of  CONT  with  a  25  x  255 
array  of  pressure  coefficients,  it  is  recommended  that  an  8  x  128  array  be 
used  with  C0NT1. 

Figure  8  is  a  wall  pressure  map  of  the  same  data  as  in  figure  7  except 
that  only  every  fourth  circumferential  point  (of  the  original  255)  was  used 
in  the  reduction  with  C0NT1  (i.e.  axial  interpolation  is  linear  through  use 
of  C0NT1  but  MAP2  is  not  used).   The  resulting  7  x  64  array  of  pressure 
coefficients  gives  rise  to  very  similar  contours  to  those  shown  in  figure  7 

for  the  finer  mesh.   In  fact,  the  appearance  of  contours  with  C  =  0.4  near 

P 

the  blades  in  figure  7  might  be  an  erroneous  result  of  quadratic  interpolation 
across  blades.  Note  that  the  grid  lines  in  figure  8  which  indicate  transducer 
locations  are  incorrectly  plotted  to  have  equal  spacing,  but  this  does  not 
affect  the  above  comments. 

Figure  9  is  a  plot  of  contours  for  a  higher  compressor  rotational  speed. 
The  data  reduction  techniques  were  identical  to  those  used  to  obtain  figure  5. 
Note  that  the  blade  location  shown  was  that  calculated  for  the  data  of  figure 
5.  The  location  is  slightly  in  error  due  to  a  difference  in  the  tuning  of  the 
phase-lock  loop  circuit  in  the  pacer  for  the  second  test.   This  problem  has 
been  solved  recently  by  using  the  near-discontinuity  in  the  pressure  distribu- 
tion measured  by  the  transducer  at  the  leading  edge  to  position  the  blade  tip. 

Again  it  is  stressed  that  the  results  presented  are  prelinimary  and  are 
intended  merely  to  demonstrate  the  methods  used  and  the  capabilities  of  the 
system. 
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5.   CONCLUSIONS  AND  RECOMMENDATIONS 

Programs  for  the  acquisition  and  reduction  of  fluctuating  casing  wall 
pressures  in  a  transonic  compressor  stage  have  been  developed  and  run 
successfully  with  the  compressor  at  this  time  operating  at  up  to  60  percent 
of  the  design  speed.   Evaluation  of  the  data  acquisition  system  on  a 
mechanical  simulator  indicates  that  it  can  operate  over  the  full  speed 
range  of  the  compressor.   In  fact,  at  higher  speeds  the  signal  to  noise 
ratio  in  the  Kulite  transducer  outputs  will  improve  significantly  because  of 
the  higher  pressures  that  will  be  encountered. 

The  pacer  mode  of  synchronized  sampling  has  made  it  possible  to  deter- 
mine in  a  versatile  manner  the  wall  pressure  maps  in  the  rotating  frame  of 
the  rotor.   The  pacer  system  can  also  be  used  to  obtain  measurements  of 
flow  properties  away  from  the  wall,  e.g.  flow  velocity  measurements  with  a 
dynamic  probe.   Wall  pressure  maps  have  been  presented  solely  to  demonstrate 
the  capabilities  of  the  pacer  technique  and  of  the  data  acquisition  and 
analysis.   Comprehensive  data  will  be  presented  and  interpreted  in  a 
subsequent  report. 

There  are  some  aspects  of  the  programs  which  can  be  refined  or  which 
need  further  evaluation  and  the  following  recommendations  are  made. 

1.  The  subroutine  RPACE  in  KULITE  causes  a  sample  to  be  taken  in  pacer  mode 
about  every  tenth  revolution  of  the  rotor  at  60  percent  of  design  speed. 
This  causes  a  delay  in  data  acquisition  which  could  be  reduced  by 
modifying  subroutine  RPACE. 

2.  The  degree  of  steadiness  of  pressure  distributions  in  the  frame  of  the 
rotor  needs  further  investigation.   This  should  begin  with  an  examination 
of  the  standard  deviations  already  computed  in  pacer  mode.   In  separate 
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tests,  a  larger  number  of  samples  (at  least  500)  should  be  taken  for  each 
of  several  steps  between  blade  pairs  and  the  variations  at  each  step 
examined  for  frequency  content.   The  measurements  should  then  be  repeated 
with  the  case  wall  rotated  peripherally  by  at  least  90°. 

3.  The  technique  of  calibration  of  the  Kulite  pressure  transducers  under 
operating  conditions  effectively  takes  account  of  change  in  transducer 
sensitivity  with  temperature.   Change  in  transducer  offset  (d.c.  level) 
with  temperature  is  presently  handled  by  equating  the  time-average  trans- 
ducer output  voltage  with  the  steady  pressure  obtained  from  a  pneumatic 
static  tap  at  the  same  axial  location.   The  relationship  between  the 
steady  pressure  indicated  by  the  pneumatic  tap  and  the  time  time-averaged 
pressure  at  the  tap  needs  further  investigation. 

4.  Because  transducer  K10  and  pneumatic  tap  S10  are  not  coincident  axially 
it  is  necessary  to  interpolate  between  readings  at  S10  and  Sll.   The 
interpolation  in  MAPI  is  at  present  linear  but  its  adequacy  has  not  been 
fully  evaluated. 

5.  The  large  pressure  gradients  in  the  axial  direction  across  rotor  blades 
are  not  resolved  well  because  of  the  limited  number  of  transducers.   Two 
additional  transducers,  located  midway  between  K9  and  K10  and  K10  and 
Kll,  would  greatly  alleviate  this  difficulty.  Linear  axial  interpolation 
across  blades,  as  in  C0NT1,  is  misleading  and  hand  smoothing  of  contours 
near  blades  is  presently  necessary.   Quadratic  interpolation,  as  in  MAP2, 
does  not  solve  the  problem.   Extrapolation  of  data  up  to  but  not  across 

a  blade  surface  should  be  investigated. 
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6.  Shock  waves  have  not  been  encountered  at  the  low  operating  speeds  at 
which  the  present  data  was  obtained.   The  accuracy  of  resolution  of 
shock  waves  should  be  studied  in  the  light  of  the  above  discussion  re- 
garding large  pressure  gradients  across  blades. 

7.  The  blades  have  been  located  on  the  wall  pressure  maps  in  this  report 
from  knowledge  of  the  point  in  the  circumferential  pressure  distribution 
(indicated  by  transducer  K9  at  the  blade  leading  edge)  at  which  the 
circumferential  pressure  gradient  is  steepest.   This  technique  is  subject 
to  an,  as  yet,  undetermined  uncertainty  due  to  irregularities  in  the 
geometry  from  blade  to  blade.   The  alternative  procedure,  whereby 
blades  are  located  by  use  of  the  phase  relationship  between  the  one  per 
blade  signal  and  the  circumferential  pressure  distributions,  also  needs 
further  evaluation. 

8.  In  its  present  form  the  data  acquisition  system  requires  frequent  key- 
board entries  by  the  operator.   In  principle  the  system  can  be  fully 
automated  by  pre-entering  all  necessary  data  with  DATA  statements  and  by 
replacing  INPUT  statements  by  READ  statements.   Some  WAIT  statements  in 
KULITE  would  be  needed  to  allow  TRAN4  to  catch  up  to  KULITE. 

9.  The  format  of  graphical  outputs  can  be  improved  by  using  the  plotter 
to  add  alphameric  information. 

10.   Two-way  data  transfer  beteen  the  HP21MX  and  the  HP9830A  is  feasible. 

This  capability  should  be  developed  to  enable  use  of  the  faster  HP21MX 
for  repetitious  data  reduction. 
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APPENDIX  A 

Table  1  contains  the  axial  and  circumferential  location  of  the  Kulite 
pressure  transducers  and  the  axial  location  of  the  pneumatic  static  pressure 
taps. 


Kulite 

Pneumatic 

Axial  distance 

Circumferential 

transducer 

static  tap 

downstream  of  K6 

location  relative 

numb  er 

numb  er 

(inches) 

to  K6  in 

direction  of  rotation 

K6 

— 

0 

0° 

K7 

— 

0.50 

+  10° 

K8 

— 

1.00 

0° 

K9 

— 

1.37 

+  10° 

K10 

— 

1.75 

0° 

Kll 

— 

2.12 

+  10° 

K12 

— 

2.50 

0° 

K13 

— 

3.00 

+  10° 

S6 

0 



S7 

0.50 



S8 

1.00 

_ 

S9 

1.37 



S10 

1.55 



Sll 

2*12 



S12 

2o50 

— 

S13 

3.00 

— 

Table  1  Location  of  Kulite  pressure  transducers  and  pneumatic  static  pressure 
taps. 
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APPENDIX  B 
DETAILS  OF  KULITE 

KULITE  is  the  data  acquisition  program  for  the  HP21MX  computer.   Its  three 
modes  of  operation  are  indicated  in  section  3.2.   Figure  10  is  a  flow  diagram 
of  the  program  and  a  listing  is  given  in  Table  2. 

Variable  Assignment  for  KULITE 

11  -  Run  #.   Same  as  Run  #  in  Log  Book  assigned  to  each 

start-up  of  the  compressor 

12  -  Test  #.   Refers  to  a  particular  operating  condition 

within  a  run. 

13  -  Day 

14  -  Month 

15  -  Year 

16  -  A/D  converter  mode  #. 

17  -  Samples /channel  in  free-run  mode. 

18  -  Not  used. 

19  -  Experiment  //.   Refers  to  either 

(i)   One  time  series  of  free-run  data, 

(ii)  A  complete  set  of  calibration  readings  (averaged) 

for  all  transducers,  or 
(iii)  Averaged  pressures  across  one  blade  pair  in  Pacer 
mode . 

Al  -  Channel  //.   Refers  to  A/D  converter. 

Tl  -  Transducer  #. 

Nl  -  Samples /point  in  Pacer  mode. 

N2  -  Blade  pair  #. 

M  -  Mean  of  pressure  samples  at  a  point  in  Pacer  mode. 

S  -  Standard  deviation  of  pressure  samples  at  a  point  in 

Pacer  mode. 

R  -  Mean  of  one/blade  signal  at  a  point  in  Pacer  mode. 

L  -  Row  number  in  K  matrix  of  calibrations. 

A3,A4,A5,A6  -  Associated  with  subroutine  RPACE  and  defined  by  WEST  (1976) 

A[101,  16]  -  Consecutive  free-run  samples  stored  row  by  row. 

B[101,  16]  -  Buffer  in  subroutine  R5610 

C[10],   D[10]  -  Buffers  in  subroutine  R5610 
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E[4,  255] 


K[5,  12] 


M, 


M_  -  -  -  -  step  2  -  -  -  M 


255 


s1   s3  s255 

A41  M3 M255 

\        R3 R255 


This  is  the  matrix  of  averaged  data  taken  in  RPACE 
across  two  blade  passages.   Note  that  128  points  are 
taken  across  two  blade  passages.   This  number  can  be 
changed  to  255  by  changing  Line  195. 

Matrix  of  calibrations.   Rows  1,  2,  3,  each  contain 
averaged  calibration  voltages  for  the  twelve  trans- 
ducers.  Each  of  rows  1,  2,  3  corresponds  to  a 
different  calibration  level.   Row  5  contains  II,  12, 
13,  14,  15,  19  and  the  three  reference  pressures Cwhich 
are  keyed  in  on  request)  in  K[5,7],  K[5,8],  K[5,9].   Row 
4  is  treated  as  another  calibration  level  and  is  used 
to  scan  the  offsets  if  needed.   In  that  case  any  value 
can  be  input  to  K[5,10]  for  P  Ref. 
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Figure  10.   Flow  diagram  for  KULITE 
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Table    2.      Listing  of  KULITE 


I  REM     "KULITE"     DATALOGGING    PROGRAM-- S IMMOMS- -  1  1     APRIL,      1977 
5      DIN    AC  1  0  1,  1*3,  BC  10  1,  16!!,CC  101,  DC  101,  EC  4,  2551, KC  5,  12] 

9  PRINT    "ENTER RUN*, MONTH, DAY, YEAR" 

10  FOR    1=1     TO    5 

II  FOR   J= 1    TO    12 

12  LET   KC  I,J]  =  0 

13  NEXT    J 

14  NEXT    I 

15  INPUT    I  I ,  I  4,  13,  I  5 

16  PRINT    "ENTER    0     IF    9830     IS    RESET 1     FOR    END" 

17  INPUT    R8 

18  IF    R8>.  5    THEN    47  0 

20       PPINT    "ENTER TEST*,     EXPERIMENT*" 

25  INPUT  12, 19 

26  PRINT  "0  FOR  CALIBRATION OTHERWISE  1" 

27  INPUT  Z  1 

28  PRINT*  BIZ  1 

29  IF  ZK.5  THEN  129 

30  PPINT  "ENTER MODE*,  CHANNEL*,  TRANSDUCER*" 

35   INPUT  I6,A1,T1 

40   PRINT*  8;  I  1 

45   PRINT*  8; 12 

50   PRINT*  8; 13 

55   PRINT*  8J 14 

60   PRINT*  8;  I  5 

65   PRINT*  8;  I  6 

70   PRINT*  8,#  19 

75   PRINT*  8,*T1 

80   PRINT*  8,*A1 

85   IF  16=0  THEN  165 

90   IF  16=4  THEN  103 

95   PRINT  "WRONG  MODE*" 

98   GOTO  3  0 

10  3       PRINT    "ENTER       SAMPLES    /    CHANNEL     (NOT    >     1616)" 

10  5       INPUT    17 

110       PRINT*    8; 17 

115  R5610C7, AC  1,  11,  17, Ai,  I  6, EC  1,  11 ) 

116  PPINT    "DATA    TAKEN" 
120       FOR    J=  1     TO     10  1 

124  FOR    1=1     TO     16 

125  PRINT*    8; AC  J,  IT 

12  6  NEXT  I 
127  NEXT  J 
123       GOTO     16 

129  LET  L=0 

130  LET  L=L+ 1 

13  1   PRINT  "ENTER  P  REFCIN.  H20  REL  ATMOS)  IF  READY  FOR  CALIBRATION' 

13  2   INPUT  KC  5, 6+L 1 
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Table    2.      Cont. 


133  FOP.  J  1=  1    TO    12 

134  PRINT    "STARTING    CALIBRATION    OF    A/D    CHANNEL    #    "J  1 

135  LET    A1=J  1 

13  6  R5610<7, AC  1,  13,  1616,A1, 4, EC  1,  13 ) 

137  LET   KCL,  J  1  3  =  0 

138  FOR   J= 1    TO    10  1 

139  FOR    1=1    TO    15    STEP    2 

140  IF    ABS<ACJ, I] )< .98    THEN     145 

141  PRINT    •**»*•**■»#•**»•*■*»**■*#*****■*■*•»•*■** " 

142  PRINT    w****»*OVERLOAD    ON    "J  1  -»»»*** - 

14  3  PRINT    "»■*•****♦■*****■*■*•*■»****■***##**•» " 

144  GOTO     13  1 

145  LET    KCL, J 1 3 =KC L, J  1 3 +  AC  J, 13 

146  NEXT    I 

147  NEXT    J 

145  LET    KCL, J  1 3=KCL, J  1 3/1616 

149  NEXT    J  1 

150  PRINT    "0    FOR    ANOTHER    CALIBRATION  —  OTHERWISE    1" 

15  1  INPUT    Z2 

152  IF    Z2<  .  5    THEN     130 

153  LET    KC  5,  13  =  1  1 

154  LET    KC  5, 23  =  12 

155  LET    KC5,  33  =  13 

156  LET    KC  5, 43  =  14 

157  LET    KC  5,  53  =  1  5 
153  LET    KC  5,  63  =  19 

159  FOR    1=1    TO    5 

160  FOR   J=l    TO     12 

161  PRINT*    8;KCI,J3 

162  NEXT   J 

163  NEXT    I 

164  GOTO     16 

165  PRINT    "ENTER    BLADE    PAIR    #,     SAMPLES/POINT" 

166  INPUT    M2,N1 

167  PRINT*    8;N1 

168  PRINT*    8;N2 
170  LET    A3=0 

180  IF    Tl>INT(Tl/2)*2+. 1     THEN     190 

18  1  LET    A6=32768+N2*256-64 

182  GOTO     195 

190  LET    A6=32768  +  N2*256 
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Table    2.      Cont. 


195  FOR    1=1     TO    255    STEP    2 

200  LET    A3  =  A6+I 

215  LET    P=0 

225  RPACECA3, A4,  A5) 

230  R5610(7,  CC  13,N1.,  Al,  0,  DC  m 

240  FOR   J=  1    TO   Nl 

241  LET    BCJ,  H  =  CCJ] 

242  NEXT    J 

243  LET    A3=32768+N2*256+I 

244  RPACECA3, A4, A5> 

245  R5610C7,  CC  13,N1,  0,  0,  DC  H) 

246  FOR   J= 1     TO    Ml 

247  LET    R=*R+CCJ3 
250  NEXT    J 

25  5  LET    R=R/N1 

260  LET    S=0 

27  0  LET    M=0 

28  0  FOR   J=  1    TO    Nl 
290  LET    M  =  m+BCJ,  1  3 
30  0  NEXT    J 

310  LET    M=M/N1 

32  0  FOR   J= 1    TO    Nl 

330  LET    S=S+(  (EC  J,  13-M)*CBCJ,  13-M)  ) 

340  NEXT    J 

350  LET    S=SQR<S/CN1- 1 )> 

360  LET    EC  1,  I]=M 

370  LET    EC2,  I  3  =  S 

380  LET    EC  3,  I3=A4 

38  5  LET    EC  4,  I3  =  R 

39  0  NEXT    I 

40  0  FOR    J= 1     TO    4 

410  FOR    1=1    TO    255    STEP    2 

420  PRINT#    8; EC  J, I 3 

43  0  NEXT    I 

44  0  NEXT    J 
450  GOTO     16 
47  0  END 
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Notes  on  KULITE 

1.   A/=D  converter  channels.   It  is  essential  that  the  "raw"  one  per  blade 

signal  be  input  to  channel  1  of  the  A/D  converter.  Allocation  of  the 

other  channels  is  not  unique  but  the  allocation  in  Table  3  is  recommended. 

Channels  11  and  12  are  scanned  but  at  present  are  not  used  in  subsequent 

analysis. 


A/D  Converter  Signal 

Channel  Number 

0  one  per  blade  raw  signal 

1  K6 

2  K7 

3  K8 

4  K9 

5  K10 

6  Kll 

7  K12 

8  K13 

* 

9  P    -  P 

ref    atmos 

10  S2  -  P   . 

ref 

11  Unused 

12  Unused 

Table  3.  Allocation  of  signals  to  A/D  converter  channels. 

P   .  is  pressure  applied  to  reference  side  of  Kulite  transducers, 
ref 

2.  Subroutine  R5610  is  described  by  WEST  (1976,  p.  17). 

3.  In  calibration  mode  the  scan  through  the  twelve  channels  must  be  made  four 
times.   The  first  scan  must  be  with.- the  pressure  tapping  S2  applied 
simultaneously  to  the  reference  side  of  Kulite  transducers.   The  second 
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and  third  scans  must  be  made  with  other  steady  pressures  applied  to  the 
reference  side  of  the  KulitesQ  The  fourth  scan  must  be  made  to  satisfy 
the  program  but  at  this  stage  the  data  taken  is  not  used  in  subsequent 
analysis.   This  scan  is  included  to  enable  logging  of  the  offsets  on  the 
Kulite  amplifiers  should  they  be  of  interest. 

4.  The  program  searches  for  overloads  (i.e.  greater  than  0.98  volts  or  less 
than  -  0.98  volts  in  the  calibration  signals).   If  it  detects  an  overload 
among  alternate  samples  in  the  1616  samples  taken  from  any  transducer  the 
offending  A/D  converter  channel  number  is  displayed,  the  scan  is  aborted 
and  the  program  is  reset  to  repeat  the  scan.   The  limit  of  0.98  volts  can 
be  changed  in  line  140.. 

5.  The  one  per  revolution  signal  from  the  timing  wheel  indicates  the  origin 
for  circumferential  measurements  around  the  rotor.   The  pacer  then  uses 
the  one  per  blade  signal  to  divide  the  rotor  circumference  into  9  equal 
intervals,  the  first  interval  beginning  at  the  origin.   These  intervals 
are  designated  by  blade  pair  numbers,  although  the  start  of  an  interval 
need  not  coincide  with  a  blade  tip.   Each  interval  represents  a  circum- 
ferential length,  in  the  rotor  frame,  equal  to  that  of  the  arc  (measured 
at  the  wall)  across  two  adjacent  blade  passages.   Each  of  the  nine 
intervals  is  divided  into  256  sub-intervals.   In  pacer  mode  the  scan 
across  an  interval  begins  after  the  first  sub-interval  and  ends  after  the 
255th  sub-interval.   With  stepping  sequentially  across  the  sub- intervals 
in  pairs,  a  total  of  128  points  are  sampled.   It  is  convenient  to  take  10 
samples  at  each  point  (one  sample  approximately  each  ten  revolutions)  to 
compute  the  mean  and  standard  deviation. 
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6.  Even  numbered  transducers  (e.g.  K6,  K.8  etc.)  are  located  on  one  axial 
line  and  odd  numbered  transducers  are  on  another  axial  line  which  is  dis- 
placed around  the  casing  wall  by  10  degrees  in  the  direction  opposite 
that  of  rotation  of  the  rotor.   The  parity  of  the  transducer  number  is 
evaluated  in  line  180. 

The  variable  A3  determines  the  time  (in  terms  of  degrees  of  rotation 
of  the  rotor)  after  the  one  per  revolution  pulse  when  a  sample  is  taken 
at  point  I.   For  example,  in  line  190,  A6  =  32768  +  N2  *  256  for  odd 
transducer  numbers,  and  A3  =  A6  +  I .   This  defines  the  sampling  time 
(approximately  each  tenth  revolution)  for  point  I  in  blade  pair  N2 . 
(I  =  1  to  255,  N2  =  1  to  9) .   Point  I  can  be  sampled  10  degrees  earlier 
for  even  numbered  transducers  by  setting 

A6  =  32768  +  N2  *  256  -  64 
and  A3  =  A6  +  I 
The  subroutine  RPACE  is  described  in  more  detail  by  WEST  (1976)  . 

7.  During  each  scan  of  a  transducer  across  a  blade  pair  the  raw  one  per  blade 
is  also  sampled.   This  is  used  later  in  TRAN4  to  determine  the  location 

of  the  measured  pressure  distribution  relative  to  the  rotor. 
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APPENDIX  C 


DETAILS  OF  TRAN4 

TRAN4  is  the  data  acquisition  program  for  the  HP  9830A  programmable  calculator. 
Figure  11  is  a  flow  diagram  and  a  listing  is  in  Table  4. 

Variable  Assignment  for  TRAN4 

-  RECORD  //.   i.e.  Number  of  first  available  record  in  DATAY1 

-  IDENTIFIER  (  =  0  FOR  CALIBRATION  -  OTHERWISE  1) 

-  BLADE  PAIR  # 

-  SAMPLES /POINT  in  Pacer  mode. 


Nl 
Zl 
A3 
A2 
Tl 


T2 
A[5,128] 


B[9] 


-  LOCATION  (between  1  and  255)  of  point  in  blade  pair  where 
1/blade  signal  is  0.5  volts  and  increasing. 

Location  of  point  in  blade  pair  where  1/blade  signal  is  0.5 
volts  and  decreasing. 


M, 


M, 


step  2 


M 


255 


255 


A4 


255 


R, 


1     3 

^      R2  ^^ 

II,  12,  13,  14,  15,  16,  17,  Tl,  Al,  A2,  A3,  0,  0,  0 


NOTE  that  variable  names  in  this  array  are  those  used  in  21MX 
program.   They  should  not  be  confused  with  variables  used  in 
TRAN4.   This  is  the  matrix  of  averaged  data  taken  in  RPACE 
across  two  blade  passages.  Note  that  128  points  can  be  changed 
by  changing  dimension  statement  and  the  FOR  loop. 

Buffer  for  identification  data 

[II,  12,  13,  14,  15,  16,  19,  Tl,  A1]T 

NOTE  that  the  variable  names  in  this  array  are  those  used  in 
21MX  program. 
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D[16,  102]    — 


Free  run  data 
Stored  row  by  row 


I  II 

I" 

'13 

I  14 

lis 

1 16 

1 19 

'  Tl 

I 

,A1 

I  17 
Column 


102 


NOTE  that  variable  names  in  column  102  are  those  used  in 
21MX  program. 

K[5,  12]   -   Same  as  K[5,  12]  in  21MX  program. 


DATA  FILES 


RECY  # 


This  file  contains  1  record.   It  contains  Nl  which  is 
number  of  the  1st  available  record  in  file  DATAY1 


the 


DATAY1 


Data  file  for  calibration,  free  run  and  Pacer  data.   ie.  for 
K,  D  and  A  arrays  respectively.   It  contains  300  records. 
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Q  Q,TAft-r) 


FftgE-  RvJt-A 


Recede-  from  2iw\x- 

t>AT£,  RuiS^    TE4T»-, 


3^oy^es 


Recede   SA^'p\-es/crt*>Nr*ei_ 


370  V 


PsSCeivS"    FREE'-P.-Jr*    SArAV>i£S 
A-t)E»    lt>efNTiPicATiOK    TO  b\j,3~l 


460r 


tM.SC. 


A.%0 


INUEM6NT    aeco^Tj   it 


90 


deceive    -z_-\    C=o 

?=oR     CAv_viJRAT<Oi-A    FRorv\  girw^ 


17 


C  A  U\e>R  AT  I  ONl 


ZOO 


RHCEi^G1    cAUraRATvOr*  t>ATA     ¥.V\X,3l 


3oo 


Store     v<.^.x,3l    om 

~t>\SC 


3,Q  Y 


iNCRG^Gi^T    R£LORb  *fc 


I 


^ 


->— ' 


Rcc^j^r       <a.l    A>sP     VA.ia.tb 


SA<v\P«_es/pOnHT  ,    PACeR 


t>ATA 


-«Orv<\     "2.\iv\x 


f  boo 


UOCATE-      l/8UAt>G-     Pvjv-S.5" 


IMSPUA^ 
"PO  USE- 
NOT     FOONtl 


At>t>    'C6>»T<PiCA"Tv0PA      "TO 
P«£R     T>A-TA     rv>A-TRi> 
Atl,3l 


79  O 


I 


Store    pacsr   b/viA   on 
t>\  sc 


>KCR6iv>eMT      RECOftb    *fc 


T 


OUTPUT     \t36im  iPlCATiOH 
Or4    VRlSTER 


SES 


Figure   11.      Flow   diagram  for   TRANA 


PACER 
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Table    4.       Listing   of   TRAN4 

i       •  • :.  [  MM  01  13 i.  2    fl!  ! 


I  "I 


■"'!':,     :::    :  i    •*  ■:.'•     p :;:"  i'-    '  >'      *  i ;.  •  ra  ,■'  i 

■   .    '.J  I         ;.   L.I..-...'  '•     :_  '...      !      :'  J  r  I     I     i    I     I       : 

40    D  [  i'1    DSC  1 ,:;:  j  1  u 2  1 1  fl  ~!t  5  i  I  28  h  t  £  9  3?  1  i  Ffi .■•  i  ':•  3 

53    REH.D    #15  HI 

60    MAT    D-2ER 

79    MAT    R---ZER 

88    REM    TEST    IF"    HE  XT    DATA    FROM    21  MX    IS    H    i.-tfj   IBPflTIul- 

9  8    E  H  f  E  R    '  i  ••  *  '.-  •: !. 

1  Ski  IF  ZK8.5  I'HEH  200 

1  1 9  R  E  M  R  E  C  E  I V  E  1 1)  E  f  I T  I F  I  C  fl  T  I  0  |  i  D  R  f  R 

120  FOR  1=1  TO  9 

1  30  EH  TEE  •  1»4  :-eC  I  1 

140  HE XT  I 

159  REN  TEST  MODE!* 

169  IF  BE  6  3=0  THEN  530 

[70  IF  BE- 6  3=4  THEN  360 

180  DISP  "WROHG  MODE  #"j 

1  i-  O  R  E  N  R  E  C  E I V  E  R  H  L  S  T  0  R  E  C  R I .  I B  R  fl  T 1 0  N  D  R  T 1 1 
3  0  0  F  0  R  I  -  1  T  0  5 

210  FOR  J=l  TO  12 
220  EH TEE  < 1 ? *>KC I  j  J  3 
230  NEXT  J 
240  NEXT  I 
250  PRINT 

2  6  0  P  R I N  T  ' '  R  R  R  R  V  0  F  C  fl  L I B  R  R  i"  1 0  N  V  0  L  T  fl  G  E  S  ' 
270  PRINT 

220  MAT  PRINT  K 

290  READ  #2? HI 

300  MAT  PRINT  #  25 K 

310  H1=N1+1 

320  READ  #1 » 1 

3 SO  PRINT  #1?N1 

3  4  0  G  0  T  0  9  2  O 

350  REM  RECEIVE  AND  STORE  FREE-RUN  DATA 

360  ENTER  (Is *)I7 

370  FOR  1=1  TO  16 

3  8  O  F  0  R  J  =  1  T  0  1 8  1 

390  ENTER  ■  1 > *)DC  E  ?  J  3 

400  NEXT  J 

410  HE XT  I 

420  FOR  1=1  TO  9 

430  JJC  I  j  102  3--BE  I  3 

440  NEXT  I 

450  DC  10 j 102  3=17 

460  READ  #2»N1 

470  HAT  PRINT  #  25  D 

430  N1=H1+13 

490  READ  #1)1 

500  PRINT  #1JN1 

510  GOTO  350 

520  REM  RECEIVE  PACER  DATA 

530  ENTER  (  1)  OA2 

540  ENTER  <1j*>A3 

550  FOR  J=l  TO  4 

560  FOR  1=1  TO  123 

5'70  ENTER  (  i  j  -..'Hi  J »  I  3 

580  NEXT  I 

590  NEXT  J 

608  REM  FIND  POSITION  OF  REFERENCE  PULSE- 1/  BLADE 
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Table   4.      Cont. 


b  1.8  i  -:6 

6  iQ  [  =  1  i-l 

638  IF    1M.  28     fHLiI    /,:■'; 

6  ••+  8  I  f    R  [  ,:-i  ••  i  J  •  8  ,.  '  ■      H  l  r  i    i    6 

6  5  0  G  0 1 0    6  2  8 

668  Tl-I 

678  1=1+1 

638  IF    RE  4p  1  V.iU  5     (i!L 'N    ?68 

690  GOTO    678 

788  T£=l 

718  60 TO    738 

?' , : 8  I)  I S P    '  RE F E F  EM C 6    F U !  S E    r ! 0  f    F  0 U r ! D  '  - 

768  FOR    1=1    TO    9 

740  RC5»  I  3=BC  I  J 

758  NEXT    I 

768  AC  5>  16  ]=A2    ' 

778  RC  5? 1 1  3= A 3 

786  REM    STORE    PACES    DAT ft 

796  READ    #2? Hi 

666  MAT    OR  I  NT    #    25  fl 

818  Nl=Nl+5 

828  READ    #1  j  L 

838  PRINT    ttlSN] 

8  46  R  E  M    P  r  i\\  T    E  X  i::  EP I  PIE  N  T    I D  E  N  T I F  A  C  T  1 0  H 
858  PRINT 

868  PRINT 

878  WR 1 TE  (  1  5 ,  686  >  BE  1  3?  BE  2  3?  BE  4  ]  j  BC  6  I s  BC  5  3 

388  FORMA  V  "X ', *  '  RUN#  "  >  F5 .  Q  j  5X  j  "  TES  T# "  <  F5  .  6  p  L  66  ■  F  3  .  0 »  "  ■■-' "  j  F3 .  6 »  " ,  "  ,  F5  .  6 

896  PR  IN T  " M 6 D 6  #  B C  6 3 > "      E X P E R I M E N T  # " B E  7 ' 3 

908  ir  BE  6  3=0  THEN  1656 

918  IF  21/6.5  THEN  996 

926  PRINT 

986  PRINT 

94  6  WR  I  TE    *:  1 5  :>  880  >  61  5  j  1  3  j  KC  b  j  2  3  s  KC  5  *  -  3  ?  KC  5  ?  3  1  p  KE  5  ?  5  3 

956  PRINT    :,EXPERIMENT#"KC5p6j 

9 66  PRI N T    "  6  A L  E  B R  A  T  I  0 H  - —  P    R E F  =  '  K E  5  j  7  J 

978  PRINT  ■  REC#"N1-1 

986  66T6  1818 

9  9  6  P  R  I N  T  ' '  T  R  A  r  i  £  D  U  i  J  E  R  # "  B  [  8  ] »  "      C  H  A  N  N  E  i. ...  # "  B  C  9  3 
1666  P R  I  NT  "  S  A M P L E S - ' 6  H  A N N EL "  1 7 »  "      R E C  tt  "  N  1  -■  1  3 
1618  PRINT 

1628  PRINT 

1646  GOTO  I  126 

1 656  PR  I  NT  '  TR  hNO  DUCER  #  ;  BE  8  3 j "      CHANNEL  # " BE  9  3  * 

1666  WR I T E  < 1 5 < 1688 ) A3 j  A3 » N 1 - 5 

1678  PRINT  "REFERENCE  PULSE  AT  POINT#"  (T  l+T2>-'2 

1 886  FORMAT  "  BLADE  PHI  R  #"  i  F4 .  8 >  5X  j  "  SRMPLES/PO I  NT  "  »  F4 .  6 ,  5X  t  "  RE6  tt "  »  F5  ,,  6 

1696  PRINT 

1168  PRINT 

1    1  l-"l      P  P  T  ri  I        '     :c  4  *+■  -lr  • '"  ***r  -+*•+■  ■+■  "*r  "A*  -it '+-'  :r  ■+■?**  ■+•  +■  r*r  ""f  '£■  H~  ?:  ;t  7-  ?:  -»r  +•  ■••■  -*r  ■*»:•  i*r  >r  -•:  ■+•  +■  :•:  'S'  '*«r  t  ■+■  x  #  r*r  "i*r  r/  ~*r  *+•  ■:•*•  ?*r  'i*r  -;tr  :*r  ir  :'*  i-  t*   )■■  '■  ■♦? 

1126  DISP   ENTER  1  F6R  MORE  1ihTA;!  j 

1138  INPUT  Q 

1146  IF  Q=l  THEN  38 

1159  PRINT 

1168  PRINT 

1170  PRINT 

H66  END 
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Notes  on  TRAN4 

1.  (line  30)  DATAY1  is  a  file  on  a  removable  disk  for  temporary  storage.   At 
the  end  of  Run  number  n   (n  is  a  two  digit  integer)  a  file  CKRWn  must 
be  opened  on  the  fixed  disk  and  data  must  be  copied  into  it  from  DAYA1 
for  long  term  storage.   The  number  of  records  in  CKRWn  must  equal  the  sum 
of  record  number  printed  out  with  identification  of  the  last  experiment 
and  a  number  k  where 

k  =  1   if  last  experiment  was  a  calibration 
=5   if  last  experiment  was  in  pacer  mode 
=  13   if  last  experiment  was  in  free- run  mode 
Figure  12  contains  sample  print-out  from  TRAN4  for  the  three  modes. 

2.  (Line  600)   The  position  of  the  centre  of  the  raw  one  per  blade  pulse  is 
found  relative  to  the  pressure  distribution  across  a  blade  pair  by 
searching  through  the  128  averaged  samples  at  each  point  to  find  the 
sample  which  first  exceeds  0.5  volts  and  the  first  subsequent  sample 
which  is  less  than  0.5  volts.   The  corresponding  point  numbers  are 
averaged.   If  the  pulse  is  not  found  (due  to  inadequate  signal  level) , 
the  program  displays  an  ERROR.   After  the  correct  one  per  blade  signal 
has  been  re-established  both  KULITE  and  TRAM  must  be  rerun. 
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Figure  12.   Sample  print-out  from  TRAN4  for  calibration,  mode  0  (pacer) 
and  mode  4  (free-run)  operation. 


RUH#  "':  !  £S  I  if  1  5         V, 

F  ■•.FERIMEH  i'tt    l 

CfU  I  BRAT  I  ON P    PEf  =  -4£.  2 

R£C#     ! 


i;    i-  +   +■  :h  +■    i.    .     .;•  -i-   i-   \-   [(■   -r-  -V.  -ij   -i-  ;!;■  >£  £  -if-  :•:.  j-  -:r  £  -i-  .+■  -if-  -h  .,.  ..,.    i-  ..,:  .,.-  +    |.    i     |.  i\i  .+   .,..  ,|.  ,:p.  .:/  .£  ^..  .+.  ...  -:,:  .£  .}:.  .(..  ..,:.  -{   ■(.  $ 


RUNtt    58     TESTtt    i  5 

NODE  tt  4  EXPERIMENT  tt  2 

TRANSDUCER  tt  o  CHRNNEL..  tt  1 

SflMF  .ES/ CHANNEL  1616  REC  #  2 


RUNtt    58      TESTtt     i  5/  24-'  1977. 

!  .  DL:  it  8  EXPERIMENT  #  3 

TRAM  DUCER  tt  6  CHRNNEL  it  I 

B I  R  D  E  P  1 1 1  R  ff  3  R  M  P  L  E  S  F;  0 1 N 1"   1 Q      R  E  C  tt    I ! 

REFERENCE  PULSE  RT  POINTtt  31 


■t  ~k  :•-  +  :■:   :~  *  ■+•  '.'  :':  '.'     r:  . :'  r.  :"  "'"  77  *  77  77 77  '?r  77  77  A-  77  77  tt  A1  A1  :<:■  7r 7.'  *  A"  77  •*:  'A'  'A1  tt  'A'  *  A'  77  r.'  w  77  'A1  77  'A"  *  W  *  *  "C  *  7:  77  -A'  tt  'A'  77  :' 


RUN#    58  TESTtt     1             5.  24-  19 

NODE  #  9  EXPERIMENT  tt  4 

TRfiuSDUCEF  tt  6  CHRNNEL  tt  1 

BLRDE  PR  IF  tt    3  SAMPLES.-  •  POINT   10      REC  tt 

REFERENCE  PULSE  RT  POINTtt  31 


.+.  .£.  .*.  ,£  .j. .+.  -i.  +. ..,-..  i_  ..,.   .*.    i:.  -,..  .+.   ,.  -r   +■  ■•{-  ■+•  -+.   ±  Jf  ■+■  •+■  -,:  •+-  ijj.  77  77  ■+•  -+•  77  77  tt  77  ■it'S  TT tt  77  A'  77  ■#■  A'  77  "A1  A'  $'  'A'  •*■  'I'?  77  tt  tt  •$•  77  77  ■*•  ~r  tt  $■  77  77  -1--  '■•  7. 


FUN  it         53  TESTtt  1 

NODE    tt    4  EXPERIMENT    tt    5 

f RAN 3 DUCER    tt    7  CHRNNEL.    tt    2 

SHMPLES. •■"CHRNNEL    1616  REC    #    25 
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APPENDIX  D 
DETAILS  OF  RESET1 

RESET1  initializes  to  1  the  number  in  the  file  named  RECY// .   This 
enables  data  acquired  at  the  start  of  a  run  to  be  stored  at  the  start  of  the 
file  named  DATAY1 .   Table  5  is  a  listing. 


Table  5.   Listing  of  RESET1 


R  E  8  c. !  J  "  -  3 1 M  M  0  N  S      S  £  ;  3    F  E  C  Y  #    T 


&  I:  I      U  ill    I        U 

30    FILES    RECV1 

4y   h  - 1 

•■■•  0  rerii  #i  j  i 

69  PRINT    ttlifi 

70  READ  #  I.  j  1 
8  0  R  E  R  D  #  1  5  fl 
90  PRINT  H 
108  END 
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APPENDIX  E 
DETAILS  OF  MAPI 

MAPI  is  used  to  compute  the  sensitivities  of  the  Kulite  pressure  trans- 
ducers and  to  convert  voltages  sampled  at  the  pressure  transducer  outputs 
to  pressure  coefficients.   Figure  13  is  a  flow  diagram  and  a  listing  is  in 
Table  6. 

Variable  Assignment  for  MAPI 

N      -   compressor  RPM 

T      -   total  temperature  (called  T  0  elsewhere)  measured  at  axial 
location  of  S2  (entered  in  degrees  F) . 

PI      -  Static  pressure  (P    )  measured  at  S2  (entered  in  inches  of 
water  absolute) . 

P2      -  Total  pressure  (P  _  )  measured  at  axial  location  of  S2 
(entered  in  inches  of  water  absolute) . 

M      -  square  of  Mach  number  at  axial  location  of  S2. 

2    2 
Ul      -  square  of  rotor  tip  speed  in  ft  /sec  . 

Al      -  square  of  speed  of  sound  at  axial  location  of  S2  in  ft/sec. 

Q      -  Reference  pressure  for  computing  pressure  coefficients. 

CI,  C2,  C3  -  First,  second  and  third  calibration  pressures  applied  to 
reference  side  of  Kulite  transducers  (inches  of  water 
relative  to  atmospheric) . 

RO      -  Record  number  for  calibration  data  on  disk. 

TO      -  Kulite  transducer  number. 

A      -  Kulite  transducer  sensitivity  in  inches  of  water  per  volt. 

PO      -  Pressure  from  pneumatic  wall  static  tap  corresponding  in  axial 
location  to  a  Kulite  transducer.   (inches  of  water  relative 
to  pressure  at  S2) . 
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Figure   13.      Flow  diagram  for  MAP   1 
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Table   6.      Listing  of   MAPI 


hi  : 


Li      >  :i 


;.;■!  0    i-: 

.  L  L3      :   :i  !":'F  ! 

40    D 
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t'    t«                 4.      1 

1 ;-  i      hi 
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f :  ,:3    D 

Q  !  i       7  i 

:";P    '  PMTFfr    R<  M 

..if    i 

100 

■ii  ij  i 

)ISP      'ENTER 

1  1  0 

INPUT     r 

129 

-  T  t-4Gti 

130 

,}.:,;.:■   -entef;  PS 

148 

HPU "i     f'l  j  FT 

150    i 

|=:<  .:  F   .-■     p  j     i   !   ,   i-1 

160    ' 

Ji=(PI*H*5- 5-   . 
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■    r   !•'  14 
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i".lj      !      I 


4  - 


'■.:>?  r:  z 
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£  :'  0    ii 0  3  U  £    'I ' 3 • 

£ 8 0    PR  INT     '8 E  N S      0 F     f '  T 8 "  =    fl "  I H .    H £ 0    ' ' 0 L  "I 

£  9  8    D  I  S  P     '  E !"  I T  E  R    \ '  U  i  "i '. ..  L    F  0  R    T "  T  8 "  I N  -  H  £  0    R  E  L    S  2 
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330    IF    ABS     f0-10'  '8.3     THEN    368 
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340    INPi.iT    P9 

350    PQ=P0  +■ '  P9-P8  ■  *8  .38  .  5  ? 

368    DISP    "ENTER    PflCER    PRESS    RECtt    FOR    f'TS? 

370    INPUT    Rl 

330  READ  #1jR1 

390  NAT  READ  #  1 i fl 

400  N  0  ■•  0 

4  i8  FOR  1  =  1  TO  J_; 

420  PC  I  3»  '  •  AC  1  t  N0*32+  I   3-CC  1  j  T8-5  1  •  *fl+P0)  ,-Q 

4:10  NEXT  I 

4  4  8  R  E  fl  H  tt  2  h  (  T  8  -  6  )  *  4  +  N  8 1  i 

450  NAT  PRINT  #  2?P 

460  LCT0-5jN0+1 ]=PC 1  J 

47  8  FOR  I=£  TO  32 

489    IF    PC  I  ]>LCT0-5jN0+1  ]    THEri    588 

498    LC T0-5»N0+1  3=PC  1  1 

588    NEXT    I 

518    UCT0-5»N0+1  ]=PC  1  J 

528    FOR    1=2    TO    32 

538    IF    PC  I  ]<IJET8-'5i.!-I8+1  3    THEN    558 

548    UCT0-5jN0+1 3=PC  1  3 

558  NEXT ■ I 

568  N0=N0+1 

578  IF  N0  -3,  5  THEN  596 
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590  T8  =  T8-+  1 

668  IF  T 0  :■  ■  1  3  .  5  T !■ i E H  6 2 8 


H  II  U 
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Table  6.   Cont. 


b  I 

I-  ;" 

6  S  8 

6  5  8 
66W 
670 

t  ?  0 

696 

7  0  9 
7  1  6 
720 
739 
740 
750 
760 
776 
7  S  0 

7  9  8 

8  8  8 

8  i  0 

8  -':  0 
8  3  0 
840 
8  5  0 


I  OP    I.    1   TO   8 

FOR  -..1-1  TO  4 

iF  LlL  3  j  J  3  U  i  HE 

U=UC  I <  J  J 

NEXT  J 

NEXT  ; 

PRIM"!  "CF  ilHX=' 

L- iEi-..Q 

FOR  [=1  TO  8 

FOR  J~l  TO  4 

I  F  LI.  i  ?  J  1  >i_  THE 
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NE>H"  J 

NEXT  I 

PRINT   OF  MlN=" 

STOP 
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END 
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Rl      -  Record  number  for  pacer  data0 

NO      -  NO  +  1  is  pacer  data  quarter  number.   The  128  samples  are 
divided  into  4  sets  of  32. 

U  -  Maximum  pressure  coefficient  in  set  across  two  blade  passages. 

L  -  Minimum  pressure  coefficient  in  set  across  two  blade  passages. 

AS[5,  128]  -  Same  as  A[5,  128]  in  TRAN4. 

C[5,  12]  -  Same  as  K[5,  12]  in  TRAN4. 

P[32]       -  Array  of  pressure  coefficients  for  one  transducer  and  in  one 
quarter  of  pressure  distribution  across  two  blade  passages. 

U[8,  4]     -  U[I,  J]  is  local  maximum  pressure  coefficient  for  transducer 
number  I  and  quarter  J. 

L[8,  4]     -  L[I,  J]  is  local  minimum  pressure  coefficient  for  transducer 
number  I  and  quarter  J. 

DATAY1      -   Same  data  file  as  in  TRAN4. 

PRESS  -  contains  pressure  coefficients.  The  32  records  contain  in 
order  quarters  1,  2,  3,  4  for  transducer  6,  quarters  1,  2, 
3,  4  for  transducer  7,  etc. 

Notes  on  MAPI 


1.  The  correct  CKRWn  file  name  for  the  run  under  consideration  must  be 
entered. 

2.  Storage  of  pressure  coefficients.   The  128  pressure  coefficients  asso- 
ciated with  each  transducer  have  been  grouped  into  4  sets  (quarters)  with 
each  set  being  stored  in  a  separate  record.   This  has  been  done  to 
facilitate  use  of  interpolation  programs  such  as  MAP2.   Interpolation 
expands  the  size  of  the  data  array  so  that  pressure  coefficients  must 

be  recalled  from  the  mass  memory  in  subsets  in  order  to  meet  the  storage 
limitations  of  the  HP9830A. 
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3.  Sensitivities  of  the  Kulite  transducers  are  computed  as  follows: 

CI,  C2,  C3  are  the  three  steady  pressures  applied  to  the  reference  side 
of  the  Kulite  transducers.   They  correspond  to  mean  output  voltages 
C[l,  TO  -  5],  C[2,  TO  -  5],  C[3,  TO  -  5]  from  Kulite  transducer  number 
TO.   The  sensitivity  of  a  transducer  (i.e.  the  slope  of  its  calibration 
curve  at  a  particular  mean  operating  temperature  and  pressure)  is  ob- 
tained from  a  least  squares  fit  of  a  straight  line  through  the  three 
points. 

4.  Calculation  of  pressure  coefficients. 

The  Kulite  transducer  output  voltages  E  are  converted  to  pressure 

coefficients  C   as  follows: 
P 

C  =  ((E  -  E)  *  A  +  P0)/Q 

P 

where  E  =  transducer  mean  output  voltage  obtained  during  calibration 
with  S2  on  reference  side  of  diaphragm. 

A  =  sensitivity  of  transducer  in  inches  of  water  per  volt. 

PO  ■  mean  wall  pressure  (at  same  axial  location  as  Kulite  trans- 
ducer) measured  with  pneumatic  tap.   (inches  of  water  relative 
to  S2). 

Q  =  reference  dynamic  pressure  (inches  of  water  absolute)  computed 
as  follows: 

The  reference  dynamic  pressure  is  expressed  in  terms  of  the  upstream 

density  p   and  the  upstream  flow  velocity  measured  in  the  rotating  frame 

of  the  rotor.  Hence 

12    2 
Q  =  |  p  (VZ  +  UZ) 

where  V  =  flow  velocity  at  station  S2  (ft/sec)  . 

U  =  rotor  tip  speed 

ttN   5.5  -.  /1N 

=  3Q  x  Yj-  ft. sec.  (1) 

N  =  rotor  RPM 
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It  follows  (noting  that  variable  names  are  not  necessarily  the  same  as 
in  the  listing  of  MAPI)  that 


1      2    2   2 
Q  =  -   YP  (MZ  +  lT/aZ) 


(2) 


where  P,  a,  M  are  static  pressure,  speed  of  sound  and  Mach  number  at 
station  S2. 


But 


a 


and 


M  = 


yRTt 

III 

-1 

2 
Y-l 

w-> 

(3) 


(4) 


where   TT  and   p   are  total  temperature  and  pressure  at  S2. 

By  introducing  (1),  (3)  and  (4)  into  (2),  Q  can  be  calculated  in  terms 

of  P,  PT  and  TT  . 

5.  Note  that  S10  and  KlO  are  not  at  the  same  axial  location.   For  purposes 
of  computing  pressure  coefficients  an  effective  mean  wall  pressure  at  KlO 
is  obtained  by  linear  interpolation  using  values  at  S10  and  Sll. 

6.  The  maximum  and  minimum  pressure  coefficients  are  computed  to  aid  in 
the  choice  of  contours  when  using  the  programs  CONT  or  C0NT1. 
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APPENDIX  F 
DETAILS  OF  MAP2 

MAP2  is  used  to  convert  the  8  x  128  array  of  measured  pressure  coeffi- 
cients to  a  29  x  128  array  through  quadratic  interpolation  in  the  axial 
direction.   A  listing  is  in  Table  7. 

Variable  Assignment  for  MAP2 

B[32]        -  temporary  storage  of  pressure  coefficients 

C[29,32]     -  array  of  interpolated  pressure  coefficients  across  one 
quarter  of  a  blade  pair. 

P[8,  32]     -  Array  of  pressure  coefficients  at  Kulite  transducer  locations 
across  one  quarter  of  blade  pair. 

X[8]         -  Axial  location  of  Kulite  transducers  downstream  of  trans- 
ducer K6.   (inches). 

Q       -  quarter  number. 

PRESS        -  Same  data  file  as  in  MAPI. 

INTER       -  File  for  storage  of  interpolated  pressure  coefficients  across 
one  quarter  of  blade  pair  (15  records). 

Notes  on  MAP2 

1.  Lagrangian  interpolation  is  used,  i.e.  if  P  ,  P   and  P  are  known  at 

x-,  x-,  x_  then 

(x-x2)  (x-xj       (x-Xj)  (x-x_) 
(x1~x2) (x1~x_)   1    (x2-x- ) (x2~x_)   2 

(x-x^)  (x-x?) 

+ —  p 

(x_-x., )  (x3-x2)   3 
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Table    7-      Listing  of  MAP2 
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2.  As  a  result  of  the  equispaced  interpolation,  the  transducer  axial  loca- 
tions will  in  general  not  coincide  with  axial  stations  in  the  interpolated 

array. 

3.  Three  interpolations  are  made  in  the  axial  direction  between  transducer 
measurements. 

4.  The  program  must  be  run  for  each  quarter  of  the  array  of  measured 
pressure  coefficients.  After  each  running  of  the  program  the  contours 
must  be  plotted  with  CONT  prior  to  running  MAP2  for  another  quarter. 

5.  Interpolation  between  X[l]  and  X[2]  is  made  with  a  quadratic  through 
pressure  coefficients  at  X[l],  X[2]  and  X[3].   Interpolation  between 
X[7]  and  X[8]  is  made  with  a  quadratic  through  X[6] ,  X[7]  and  X[8] . 
Interpolation  between  X[I]  and  X[I  +  1]  (for  I  >  1  and  <  6)  is  made  by 
averaging  the  quadratic  through  X[I  -  1],  X[I]  and  X[I  +  1]  and  the 
quadratic  through  X[I],  X[I  +  1]  and  X[  1+2]. 
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APPENDIX  G 
DETAILS  OF  CONT1 

CONT1  is  used  to  plot  contours  of  constant  casing  wall  pressure  from 
an  array  of  measured  pressure  coefficients.   The  program  handles  the  non- 
uniform axial  spacing  of  the  Kulite  transducers  but  requires  uniform 
circumferential  spacing  in  the  array.   C0NT1  is  written  to  accept  the  array 
generated  by  MAPI.   Figure  14  is  a  flow  diagram  and  a  listing  is  in  Table  8. 

Variable  Assignment  for  C0NT1 


P[I,  J]  -  array  of  pressure  coefficients 

A[3,  3],  B[3,  3],  Z[3],  Q[3],  D[2,  2],  F[2],  G[2,  2],  H[2]   -  defined  in 

note  6  of  this  appendix. 

x[9]  -  axial  location  (inches)  of  Kulite  transducers  downstream 

of  K6. 

10,  JO  -  dimensions  of  P[I,  J]  in  axial  and  circumferential  directic 

respectively. 

C  -  value  of  pressure  coefficient  on  a  contour. 

E  -  triangular  element  number. 

El  -  number  of  starting  element  in  a  contour  plot. 

E2  -  number  of  finishing  element  in  a  contour  plot. 

PRESS  -  Same  file  as  in  MAPI. 

Notes  on  C0NT1 

1.  The  X  (axial)  and  Y  (circumferential)  dimensions  of  P  must  be  entered 
before  program  is  run.   P[8,  128]  is  nearly  the  maximum  array  size  that 
can  be  stored  in  the  HP9830A. 

2.  The  axes  are  drawn  so  that  x  (axial)  runs  from  K6  to  K13  and  Y  (circum- 
ferential) runs  from  the  start  to  the  end  of  a  blade  pair.   Contours 
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Figure  14.   Flow  diagram  for  C0NT1 
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Table    8.      Listing  of  C0NT1 
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160  INPUT  I0jJ0 
1 7  0  S  C  A  L  E  9  n  X  L  8  3  j  9 .  5  ?  J  9  +  9 .  5 

150  XflXIS  9.5j4j0jXE83 
1 9  0  X  fi  X  1 3  J  0 +-  9 .  5  ■  3  ?  9  j  X  E  8  J 
299  • '  fl X I S  X E 1 3  n 1 3  0  j  0 . 5  j J  9  +  9 
219  V  fl  X 1 8  X  C  8  3 s  1  3  9  j  0 ,.  5  ?  J  0  +  9 

229  FOR  1=2  TO  7 

2 3 0  V fl !::l  I S  XE  I  3  ?  1 3 0  j  9 .  5  j  J  9  +■  9 
240  PEN 

259  NEXT  I 

260  READ  #1 1  1 
270  MAT  READ  #  15P 
230  DISP   ENTER  CONTOUR  P=" 
299  INPUT  C 

299  DISP  "ENTER  FIRST  & 
316  INPUT  E1jE2 
329  FOR  E  =  E1  TO  2*< 19-1 
339  IF    E/E2  THEN  1250 
349  IF  E/INT<E/2.:i*2  THEN  430 
3 59  AC  1 j 1 3  = I N T ( ( E / 2-0.90 0 0 1  > 
369  AC  1>'23=E  2- (AC  it  13- .1  >*<Jt 
379  AC  2 j 1  J -HE  1  "  1  1+1 

339  HE  2.1  2  J™  AC  1  j  2  3-1 
399  AC3j  1  3=AC  1?  1  3+1 
409  AE  3)  2  3= AC  1  t  2  3 

420  GOTO    500 

4  30  HE  1  j  1  3=  I  H  r '  ■  ■  E+  1  >  ''2-9 .  0900  1  )    ■ 

449  AC  1 .-  2  3=  (  E  +  1  ;  .-'2-  (.  AC  1  ?  J  3-  1  )  *  ;    J9 

459  AC  2 j  1  3=flC 1?  L 3 

4  60  HE  2?  2  J -HE  1.  ,  2  ]+1. 


-Mo  i     t  Lc.nfc.ri  i  w 
S-  i  ._li/|~1  '.i 


ll-i  -  I 

'I   +-   i 


-hi 


4  79 

AE3j  1  ]=flC  J  j 

1  3+1 

489 

AE  3  j  2  3=A(  !.  n 

"i  'i 

490 

R=i 

500 

AC  1^3  3=1 

5  :  O 

HE  2s  3  3=1 

520 

HE 3j 3 3=1 

5  :!0 

ZE  1.  3=PC  AE  1  • 

I    !  '■  H  L   L     ,  -  .i 

l 

•J  t  tJ 

ZC2  3=PC  AC  2 

1  i ::  r :  i   "  ;  2  ] 

J 

550 

TCP  ]  =  P[  fl[ 

I    h  Al    :  .2  ! 

! 

560 

REM    CHECI 

F     P     i  S  "IN 

!" 

570 

[F    ZE  .1  J 

L-     Al  i  iJ     1 

....'     . ' 

j    7  . ; ! 

500 

IF    ZE2  3 

I    A  N  D    L 

]     FH 

5  90 

IF    ZE  !    !      = 

r.  hnd  c 

■  zr  ' 

"i   r:-iEN 

690 

[  1     z\       3 

0    AND    C 

'  !"     ' 

I     TNI:     : 
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Table    8.      Cont. 


6  i  0    G  0  T  iD    1  2  3  Q 

629  MAT    B=INV<R) 

630  MAT    Q  =  B*Z 

640    IF    R<0.5    THEN    940 
650    IF    QC 1  ]  =  0    THEN    720 

6  6  0    X  =  '•  -  Q  C  2  ]  *  A  C  1  s  2  3  -  Q  [  3  ]  +  C  )  /  Q  C  1  3 
670    Y=AE  ls23 

680    IF    fit  1  j  1  ]    <=    X    AND    RC3sl3    >=    X 

690    GOTO    720 

700    W=INT<X> 

?  10    PLOT    X [  W  1  +  ( X [  W  + 1  3- X L  W  :i )  *  <  X - N )  j  Y 

720    IF    QC2  3=9    THEN    790 

730    Y=£-QE  1  3#RC  1 s  i  3-Qr  3  3+C)/QC  2  1 

740    X=flE  Is  1  3 

750    IF    ACls23    <  =    V    AND    AC  2s  2  3    >=    Y 

760    GOTO    790 

770    W=INT' X) 

7  s 0  p  l  o t   xc  w  i  f  < :•: C  W  +  I  ] - X C  U  3 )  *  •  X  -U ) )  Y 
790    DC  1 5  1  3=1 

800    DC  1 j 2 3=1 

818    DC  2s  1  ]=QC  1  1 

820    DC2s2  ]=QC2  3 

330    FC  1 3=RC 3s i  3+flC 3s 2  1 

340    FC2  3=C-QC  3  1 

350    n=DETCH) 

368  [F  M=Q  THEN  1238 

878  MAT  G=IHV<DJ 

888  NAT  H=G*F 

898    IF    HC 1 3    <=    AC  3s 1  1    AND    AC  1  1     >-    A 

908    GO TO    i 230 

910    14=  INK  HE 1 3) 


HEN 


f'-M 


;  L  N  +  i  j 


;  r  !,j  i 


M.  i 


920    PLOT    XCW3  + 

930    GOTO    1230 

940    [I-    QC  1  3=8    ["HEM    1018 

950    X=  (.  -QC  2  3* AC  1  s  2  3 -QC  3  3  n; )  •■•  QC  i  ] 

960    Y-AC 1 s  2  3 

970     [F    AC  1*  1  3    ■  =    X    AND    AC  3s  1  1    >  =    ! 

988    GOTO     L810 

990  W=IN  !''•;:  i 

1000    PLO r    XC N  3+  . XC N  *  1 J-XC N  3) s'X-U 

1010    IF    QC 2 3=8    THEN    1038 


1020  Y=(  -Q[  1  3* Hi.  , 

1830  X=AC 2s 1  3 

1040  [F    AC£s2  3 

1850  G 0  T 0    1 0 8 8 

1068  l-i~  I  NT1:  X 

1070  PLO  l     ':'[.  W  ]  I-1..; 

1088  DC  [  ?  I.  3=  I. 

1  8  9  0  D  C  I  s  2  1  =  I 

1  !  00  DC  2  s  1  ]=Q[  j   I 

1 i 10  DC 2i2 3=QC 2  J 

1  i  20  FT  1  ]=AC  2  s  I.  3 

1  i.  38  FC  2  3  =  C-QC  J  ! 

1  148  M=DE f  .  D) 

1150  IF    M=0    THEM 


i  J  -QC  3  3  +  C )    Q! 
'•!■'    AND    AC  3s  2  3 

'•4+1  3-Xi  W  3  '  -    ' 


.>  i  j 


3    T! 


i  i-l  r.  rH 


160    I'll 


1  i.  r8 
1  1  88 


r  r- 

i.  r 


Li  ::::  L 

H=G 

ic  i  : 

!      i  ,: 


M   ':,'     > 


80     A™  Ml  f     HE    ! 
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Table    8.      Cont 


i 
1 
1 

i 
i. 

1-2 

■i 
i 

1 

1 


1 6    PLOT    XE  W  ]+  ( XE  W+  1  ]-XC  N  ]  >  *  ( HE  1  ]-  U  .-■  j  HE  2  ] 
29    GOTO    1230 


5 
8f 


3  PEN 

i  NEXT  E 

i  DISP  "ANOTHER  CONTOUR?  1-YESj  O 

i  INPUT  C8 

IF  C8>0.5  THEN  230 
i  END 
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will  not  quite  go  to  these  extremities  because  the  pacer  system  starts 
sampling  at  1/256  th  of  a  blade  pair  and  finishes  sampling  at  255/256th 
of  a  blade  pair.   Recall  that  the  start  of  a  blade  pair  need  not  coin- 
cide with  a  blade  because  of  phase  lags  in  the  pacer  system  and  the 
location  of  the  timing  disk  relative  to  the  rotor.   Grid  lines  parallel 
to  the  Y-axis  are  the  lines  scanned  by  the  Kulite  transducers. 

3.  If  other  transducer  locations  are  used  their  axial  distance  downstream 
of  K6  must  be  entered  before  the  program  is  run. 

4.  The  variable  plotted  in  the  Y-direction  is  the  column  number  in  P.   The 
variable  plotted  in  the  X-direction  is  the  axial  location  of  a  transducer 
(downstream  of  K6)  and  is  derived  from  the  corresponding  row  number  in  P. 
When  setting  up  the  X-Y  plotter  it  is  advisable  to  make  both  the  X  and  Y 
scales  equal  to  twice  full  scale.  Note  that  the  circumferential  distance 
along  the  wall  across  two  blade  passages  is  3.847  inches. 

5.  Entering  of  the  first  and  last  element  numbers  enables  faster  plotting 
of  a  contour  which  is  known  in  advance  to  cover  only  a  limited  part  of 
the  field. 

6.  Triangular  element  representation  of  the  surface  defining  the  pressure 
coefficient  distribution,  Cp  (X,Y) 

The  surface  C  =  C   (X,Y)  is  approximated  by  triangular  elements  as 
P    P 

illustrated  in  figure  15.  Element  numbers  E  are  as  shown.   This  process 

represents  linear  interpolation  between  measured  pressure  coefficients. 

Contours  are  obtained  from  the  intersection  of  planes  C  =  constant  with 

P 

this  approximation  to  the  pressure  distribution.,   Thus  the  contours  are 
composed  of  straight  line  segments.  The  nodes  correspond  in  the  Y-direction 
(circumferential)  to  points  at  which  pacer  data  is  available.   The  node 
numbering,  in  the  local  sense,  for  typical  odd  and  even  numbered  elements 
is  shown  in  figure  15. 
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c 

o 

X- 
o 

Jr 
o 

0 

3 


c 
I 

V 
0 


x-l  =  Ati,Q 


X=3  =  AU,  21 


X  =  1  »  AL2,Q 


^ALS,ll-(J-At3,2l) 


=•  AL2,ll-(T-AL2,2l^ 


Figure  15  Numbering  of  triangular  elements  in  CONT  and  C0NT1.   Also 

numbering  of  nodes  (vertices  of  triangles)  for  odd  and  even 
numbered  elements. 
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XI 

Yl 

1 

a 

X2 

Y2 

1 

b 

X3 

Y3 

1_ 

.c_ 

If  Z  -  aX  +  bY  +  c  is  the  plane  containing  the  three  nodes  of  a 
triangle  then 

Z(l) 
Z(2) 
Z(3U 

where  Z(I)  is  the  pressure  coefficient  at  local  node  number  I  (i.e.  at 
XI,  YI). 

This  matrix    equation  is  written  as 

Z  =  Aa 
and  hence 

Q  =  A~  Z  =  BZ    say. 

If  E  is  odd,  the  node  coordinates  are  related  to  the  element  number 
as  follows: 

XI  =  A[l,  1]  =  INT  (((E  +  l)/2  -  .00001)/(J0  -  1) )  +  1 
Yl  =  (E  +  l)/2  -  (A[l,  1]  -1)  *  (JO  -1) 
X2  =  XI   ,   Y2  =  Yl  +  1 
X3  -  XI  +  1   ,   Y3  =  Yl 
Note  that  the  number  .00001  is  included  to  avoid  problems  associated  with 
round-off  error. 

Z  =  aX  +  bY  +  c  is  plane  containing  triangle.   For  Z  =  c',  line  of 
intersection  with  triangle  (i.e.  contour)  is  given  by  aX  +  bY  +  (c  -c ')  =  0, 
On  the  side  X  =  I  =  A[l,  1],    (See  Figure  15) 
Y  =  (-  aA[l,  1]  -  (c  -  c'))/b 

=  (-  Q[l]  *  A[l,  1]  -  Q(3)  +  c')/Q(2) 
provided  Q(2)  ^  0 
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If  Q(2)  =  0  ,  the  contour  is  parallel  to  the  side  in  question.   On  the 

side  Y  =  J  =  A[l,  2] 

X  =  (-Q[2]  *  A[l,  2]  -  Q(3)  +  c')/Q(l) 
provided  Q(l)  £   0 


On  the 

hyp< 

Jtenuse 

X  +  Y  =  A[3,  1]  +  A[3,  2] 

aX 

+  bY  =  c'  -  c 

Hence 

1 

1 

X 

m     A[3,  1]  +  A[3, 

2] 

a 

b. 

x 

I     c1  -  c 

■ 

This  matrix  equation  is  written  as 


-1 


DH  =  F 
Hence        H  =  GF   where   G  =  D 

If  E  is  even,  the  node  coordinates  are  related  to  the  element  number 

as  follows: 

XI  =  A[l,  1]  =  INT  ((E/2  -  .00001) /(JO  -  1))  +  1 

Yl  =  A[l,  2]  =  E/2  -  (A[l,  1]  -  1)  *  (JO  -  1)  +  1 

X2  =  XI  +  1   ,   Y2  =  Yl  -  1 

X3  =  XI  +  1   ,   Y3  =  Yl 

Z  =  aX  +  bY  +  c  is  plane  containing  the  triangle.   For  Z  =  c'  ,  line  of 

intersection  (i.e.  contour)  is  given  by  aX  +  bY  +  (c  -  c')  =  0  . 

On  the  side  X  =  I  =  A[2,  1]  , 

Y  =  (-  Q(l)  *  A[2,  1]  -  Q(3)  +  c')/Q(2) 
provided  Q(2)  t  0 

On  the  side  Y  =  J  =  A[l,  2] 

X  =  (-  Q(2)  *  A  [1,  2]  -  Q(3)  +  c')/Q(l) 

provided  Q(l)  ±   0  . 
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On  the  hypotenuse 

X  +  Y  =  A[2,  1]  +  A[2,  2] 

aX  +  bY  =  c'  -  c 
Hence 


1   1 
a   b 


A[2,  1]  +  A[2,  2] 


c  -  c 


This  matrix  equation  is  written  as 
DH  =  F 


Hence 


H  =  GF  where  G  =  D 


-1 


The  intersections  of  the  contour  with  the  sides  of  the  triangular 
elements  are  computed  in  the  manner  outlined  above.   Tests  are  performed 
to  identify  those  points  of  intersection  which  lie  on  side  of  the  triangle 
as  opposed  to  intersections  which  lie  on  extrapolations  of  the  sides. 
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APPENDIX  H 


DETAILS  OF  CONT 


CONT  is  used  to  plot  contours  of  constant  casing  wall  pressure  from 
an  array  of  pressure  coefficients  in  which  both  the  axial  and  the  circum- 
ferential spacings  are  uniform.   It  can  be  used  to  plot  contours  for  an 
array  output  by  MAP2.   The  program  differs  from  C0NT1  only  in  its  plotting 
of  x-coordinates  of  contours.  A  listing  is  in  Table  9. 

Notes  on  CONT 

1.  If  MAP2  is  used  for  interpolation,  the  file  PRESS  must  be  replaced  by 
file  INTER  and  P  must  have  dimensions  P[29,  32]  .   Contours  must  be 
plotted  by  processing  the  array  of  interpolated  pressure  coefficients 
in  four  sections.   This  requires  shifting  of  origin  on  the  x-y  plotter 
for  each  quarter. 

2.  The  variables  plotted  in  the  x  and  y  directions  are  respectively 

the  row  number  and  column  number  of  the  sub-array  of  pressure  coefficients. 
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Table   9.      Listing  of  CONT 


M    Wfll 


PRE SSL 


MAP    PROGRAM 


!  j ,      , 


15  UNIT    0 

16  FILES    INTER 

20  PEN  #39  DIM  P  105 J0) 

38  D I  M  PC  25  s  32  ]  >  PC  3  s  3  3  s  BE  3 » 3  ]  j  ZC  3  ]  s  QC  "' 

48  DISP  "ENTER  X  AND  Y  DIMENSIONS"? 

50  INPUT  10 i J0 

b  0  S C ALE  1 s 1 0  j 0  ? J 0 

70  XRXIS  0?  c;  10-1  ).-'-105  l,  10 

s 0  x a x is  j 0 f  ( 1 0 - 1  >■••'■  1 0 j  1 1 1 0 

90  VAX IS  1  j  100"  Oh J0 

108  YRXIS  10"  108?  0:>  JO 

118  READ  #1j 1 

120  MAT  READ  #  1 5 P 

170  DISP  "ENTER  CON TOUR  P=MS 

138  INPUT  C 

185  DISP  "ENTER  FIRST  8:  LAST  ELEMENT*- 

186  INPUT  ElsE2 

198  FOR    E = E 1    T  0    2  *  ( I  @  - 1  >  *  ( - 1 0 1 ) 

195  IF    E>E2    THEN    1606 

288  CF    E>IHT<E-'2>*2    THEN    298 

2 10  AC  1 5  1  3  =  I  N  T  (  (  E    2-0  „  0000  1  '     (  JO  ■  i  :     >-  i. 

220  AC  1  s  2  3=E    2-  (  AC  i  s  1  ]- 1)  *  (  J8--  1  >  + 1 

230  AC  2' 1 3=AC  1 s  1  1+1 

240  AC  2"  2  ]  =  fl[  1  j£  I ■- 1 

250  AC  3' 1  ]=fl[ 1 j 1  ]+l 

260  AC  3- 2  3  =  AC  1 t 2  J 

270  R~^ 

230  GOTO    360 

290  AC  1  j  1  3=  I  NT  t  <    E  f  1  .'  ■  2-0  .  8888  1  > 

388  AC  1j2]=  (  E+  1  :     2 - '.  AC  1  ?  1  ]-■  1  ;  *  •    :  ! 

318  AC  2?  i  3=AC  1  j  1  J 

328  AC  2*  2  3=AC  1 >2  1+1 

338  AC  3 ? 1 3=AC  Is  1  3+1 

348  AC  3s  2  3= AC  1  s  2  1 

358  R=l 

368  AC  1 j 3 3=1 

378  AC  2s 3 3  =  1 

388  AC  3s 3  1=1      . 

488  ZC  1  ]~P[  AC  1  s  i.  Is  AC  1  >  2  ]  ] 

418  ZC  2]=PC  AC  2  s  ]  3s  AC  2  s  2  J  J 

420  ZC  3  ]=PC  AC  3  j  1   is  AC  3s  2  3  I 
[F    C    IS    IN    E 


421  REN  CHECJ 

422  IF  ZC 1 3 


424  L'F 

425  EF 


j  j 
. ,  ] 

-  23  GOTO  !  30'" 

:,:  -  NAT  :   IN'- 

430  NAT  G  =  B*2 

448  :-  R  0.5 

445  IF    0i  !.  J-8 

458  :-<-•'  ■•;■![  2  3-r 

460  '-AC  I  ''  -'  3 

470  L!:"'  AC  3  *  1  J 

488  ■....OTO  588 

i-98  !   IT  Ms  Y--  i 

•  i  M-i  ■  r  ni  '■  i  ™i.'-! 


C  AND  C 

:.  AND  C 

f  1-iMD  C 

AND  C 


7i"  : 


:r 


r 


i  h ;  i 


I  it:  I 


!  -I  !   ,  i  i 


■E.:  EMBER 


2  3?HC2 


H  ;i 


a- 


Table    9.      Cont, 


AND    RC 3 ? 1 j 


5y5  'i*=  ..;-QC  1  3*AE  1  t  1  3-QC  3  3+O/QC  2  ] 

510  ft=R[  1*13 

528  IF    RCli-23    <=    V    AND    AE2<2]    >= 

530  GOTO    578 

540  PLOT    XpY-1 

570  DC  It  1  1=1 

580  DC  1  »23=1 

590  DC2p  1  3=QC  1  ]" 

bOO  BC2<2  3=QC2  3 

6i0    FC  1  3=AC3>  1  3+RC3»2  3 
620    FC.2  3=C~QC3  3 

625  M=DET<D) 

626  IF  N=0  THEN  1000 
630  NAT  OIH'/aO 

640  MAT  H=G*F 

645  IF    HC  1  3    <=    AC  3?  1  3    AND    HE  1  3    >= 

6 4 6  G 0 TO    1OO0 

650    PLOT    HC  1  3»HC  2  3-1 

660    GOTO    1O0O 

670    IF    QC 1 3=0    THEN    750 

690    Y, ~ < -QC23*AC  1-23  - Q C  3  3  +  C ) ••"' Q C  1  3 

700    Y=AC 1 ? 2  3 

710    IF    RC If  1 3    <= 

720    GOTO    750 

730    PLOT    XpY-1 

7  5  O    I F    Q  C  2  3 = 0    T  H  E  N    8  3  O 

760    Y=<-QC  1  J*AC2j  1 3-QC 3 3+C 

770    ;::;=RC  3-13 

7S0    IF    AC2<2  3    <=    ' 

7  9  O    G  0  T  0    8  3  O 

800    PLOT    ft?  Y-t 

330    DC  1 i 1 3=1 

840    D[  1,2  3=1 

850    DC  2 j 1 3=QC 1  3 

360    DC2j2  3=QC2  3 

878    FC  1  3=AC2>  1  ]+flt 

838    FC2  3=C-QC3  3 

885    M=DET<D) 

836  IF  M=9  THEN  1008 

890  MOT  G=INV<D> 

900  MAT  H=G*F  ■ 

905  IF  HC 1 3  >=  RC 1? 1 3  HMD  HC 1 

906  GOTO  1 O0O 
910  PLOT  HC  1  jjHC23-:L 
920  GOTO  100O 
1 0  08  P  E N 
1005  NEXT  £ 

100b  DIGP  "ANOTHER  CONTOUR?  1 
1087  INPUT  Co 
1 888  IF  C8 >8« 5  THEN  1 78 
1810  END 


i  t.  N 


40 


iC2 


HF  '    6 


AND  Hi 


1 


HF.N 


R[ 
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APPENDIX  I 

DETAILS  OF  PLOTS  A  AND  PLOTS  B 

PLOTS  A  and  PLOTS  B  are  used  to  plot  pacer  raw  data  and  free-run 
raw  data  respectively  against  circumferential  distance.   The  programs  are 
listed  in  Tables  10  and  11.  Note  that  PLOTS  A,  with  I  =  1  in  line  110, 
plots  the  output  of  the  designated  Kulite  transducer.   If  I  =  4  the  raw 
one  per  blade  signal  on  A/D  converter  channel  number  0  is  plotted. 
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Table  10.  Listing  of  PLOTSA 


1  ki  R  E  M *  *  *  *  *  P  L  0  T  S  R  *  *  *  *  * F: .  p  .  S  H  P  t. 

20  UNIT  U 

33  FILES  RECVtt? BATRY1 

49  DIM  DSL  5 j 128] 

50.  DISP  "ENTER  REC#" 5 

60  INPUT  H 

73  READ  #£jH 

80  MAT  RERD  ft  25  D 

90  SCALE  0  j 300 s - 1 , 5  ? 1.5 

10O  DISP  '"DRAW  AXES?   1=VES   0=NO"  5 

118  INPUT  N0 

120  IF  N0-Q  THEN  150 

1  30  XPN  IS  0 f 50j  0? 300 

148  VAN  1 3  0?  0„ 5?  -1.5?  1.5 

150  1=1 

160  FOR  J=l  TO  123 

L70  PLOT  2*, J" DC  I  -  J  1 

130  NEXT  J 

190  PEN 

230  END 
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Table  11.   Listing  of  PLOTSB 


5  R  EM   "PL  0  T  S  B "  -  S I M  M  0 N  S  -  - 2  8  J A N U  A  F:  V 

10  UNIT  @ 

20  FILES  RECY#>  DRTflYi 

38  DIM  DSC  16?  i02 J 

46  DISP  "ENTER  REC#" 5 

58  INPUT  N 

68  READ  #£?H 

78  MRT  REflD  #  21 D 

80  SCALE  8 j 17001 t1.5»  1.  5 

83  DISP  "DRAW  AXES?   1=YES   8=N0M « 

82  INPUT  N0 

83  IF  N6=e  THEN  110 

3 0  X A !::! IS  0 j  1 0 9 ?  8 j  1 7 8 8 

108  VAX  IS  0j8,5j-1.5j1.5 

118  FOR  1=1  TO  16 

128  FOR  J=l  TO  101 

1  3  8  f  L  0  T  (.1-1  >  *  1 8  l  >  J  -  D  l  I  j  J  ] 

140  NEXT  J 

158  NEXT  I 

155  PEN 

160  END 


69 


APPENDIX  J 
DETAILS  OF  TITIPK 

TITIPK  is  used  to  superimpose  blade  tip  profiles  on  wall  pressure 
contours.   The  program  is  loaded  onto  the  programmable  keys  of  the  HP9830A. 
The  program  is  listed  in  Table  12.   The  blade  tip  profile  is  tabulated  by 
PAIGE  (1976) ,  figure  2  and  table  1) .   The  axial  and  circumferential  units 
are  inches  and  the  program  is  compatible  with  a  wall  pressure  map  which  has 
dimensions  of  3  inches  axially  by  3.847  inches  circumferentially .   The 
"lower  left"  and  "upper  right"  on  the  plotter  should  be  set  to  the  corre- 
sponding points  on  the  wall  pressure  map. 

The  key  programs  should  be  "continued"  after  C0NT1  has  been  run  so 
that  the  Kulite  data  is  available  in  main  memory.   When  <CONT>  <f - >  is 
issued,  the  location  (YO)  of  the  blade  leading  edge  from  the  lower  boundary 
of  the  wall  maps  is  calculated  from  the  data  of  transducer  K9  and  appears  in 
the  display.   <CONT>  <f _>  is  then  issued  and  YO  is  requested  as  an  input. 
The  blade  profiles  are  then  drawn  and  the  key  program  ends.   The  contour 
plotting  with  C0NT1  can  be  continued  by  issuing  <CONT>  1250. 
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Table  12.   Listing  of  T1T1PK 

-■  ■ "  ! .  H:  '■■  •■■ ■ :'; -  -f .  r.  ■  >\  if :  !  ■  ■  ■  '.  i  ■  "? 

pf  03  ;:HM     :  is    I  [  •    i     t-:,  fll'f    :-:  h  Mi    !  j  fc    f,,;    NP!  |     PRESSURE 


33    r-l  i  -0 

40    3CHLE    iJ>   ;  ■  0  ■  J,  3-1  .:' 

50  DR'I  fi  ..'  .  633  ?  1  6 ..  7  ■■  >  6 „  885  i  59  ,  826 ,:!  •  3  .  3 

60  RERD  l  0  j  R0  :i  R  i  :•  G  j  ■  !  !6 

70    &  [SP    '  ENTER    Ytf"J 

80    INPUT    Y9 

90    Y8-Y9 

100    DEG 

1  10    fil  =  U  n    2~R1  )     '  R8-R1  ':■ 

1  20    fll=flTH«  fll    SQF:  :  1  -Hi  I' 2,  > 

130  Dl=fll+&1 -90 

146  X3-~8 

150  Y3-R1 

1  g  g  ': ;  4  =  ;:.;; ;;,  +  ■:;  £  g  -  3  *  p  ].  .  T  c  g  3  £  \ 

17y  Y4=Y3+<  C8~2*R1  )*SING1 

180  X  5  -  '.■  ■[  3  +  ■•  R  0  -  R  1 '  •  *  C  0  S  D  1 

190  Y5=Y3+<R8-R1 :<*SIND1 

L  0  0  R  E  M ■  ■-  •  -  P  F:  3  FILE  T  I P  <  fl  B  ''  -  - 

210  T2=-G1 

230  T 3= 90- III 

233  r4=<T3--T2)  ££1 

240  FOR  71  =  72  TO  73  STEP  ^"4 

250  1  =X  3'  F  1  *S  I  NT  1 

260  Y=Y 3--R I *C0ST 1 

2  7  0  [  f  '  r  <•■ '  •  0  •  3 .  3  4  7  f  HEN  3 1 S 
280  IF  Y+Y0<0  THEN  310 
2:-'u  PLOT  X  +  ;  !0  ?  Y+YQ 

339  GO  fO  :::20 
313  PEN 
32U  NEXT  T 1 

3  3  O  R  EN-  - S  Li  C  T 1 0  i !  3  I  D  E  ■  B  C ) 

343  72  =  -lJJ 

353  f3=2*Rl-Bl 

3  5  g  y  4  =  1  r  j  -  r  2 ;  2  8  0 

370  FOR  71-72  TO  13  STEP  14 

360  X=X5  -R0*COST  1 

390  Y=Y5+R0*SINT  I 

400  [  F  Y  +  V 3  >  3  .  3 4  7  I"  HE  N  440 

4i0  [F  Y+Y0  0  THEN  440 

420  PLOT  X  +  XS? Y+Y0 

430  GOTO  450 

440  PEN 

459  NEXT  :  I 

460  REM PROF  I LE  I.E. (CD  > 

470  T2=-bl+fll 

480  T3=i30-Gl 

4  90  T4=n'3-T2)  28 

5O0  FOR  T1=T2  TO  T'3  STEP  T4 
510  X=X4-fRl*3INTi 
570  Y=Y4+R L*C0ST1 

5  3  O  1  F  Y  +•  Y  8 )  3  .  8  4  7'  T  H  E  H  5  7  O 
543  IF  Y+Y0<O  THEN  570 
558  PLOT  X  +  X0 j Y+Y8 

568  GOTO  530 
578  PEN 
580  NEXT  n 
590  X1=X 
688  Vl=Y 
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Table    12.      Cont. 


t    L  U 

i-:  28 

6  3  8 

6  5  8 
6  6  9 
67Q 

6  8  8 
6 9  0 

7  0  0 

718 
7'  .1 8 
7  78 
740 
"i  58 
760 
770 
788 


REM- PR  EG  ,;n-  .■  SI  DE  :i»M 

•;--,-■  :  -Rl  -  [HlM 

f  or  i  i  n:  20 


i  ""  L  ,  :- 


■■...,  -:  ;:ou 


IF  Y+Y8.>2»  847  THEN  7 

IF  Y+Y0<@  ["HEN  "'iC 

PLOT  N  +  xOn  •;  +V0 

GOTO  720 

7' EN 

NEXT  I 

PEN 

IF  Hi) 1.5  THEN  "88 

N1--H1  +  1 

■f'0=Y9+l  ,  32  35  .-■•.  7-2  *N1 

GOTO  280 

END 


i  8  REf'i  -- --ROU  V  [HE 

30  W=3.347 

48  FOR  I  =  i.  TO  128 

50  IF  PC  4j  I  K8  THEN 

60  V1=PC 4?  I  J 

70  GOTO  -.40 

SO  IF  1=1  THEM  148 

98  IF  VK0  THEN  148 

100  V2=PC 4^  I  1 

1 18  Y0=W*<  I-i  .  5-r-Yl  ■ 

1  20  7i ISP  " Y8-  ' YO 

1  38  STOP 

'  iO  NEXT  f 

158  END 


FIND  BLADE  LEADING  EDGE  FOR  J 


1-V2 
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APPENDIX  K 

OPERATING  PROCEDURE  FOR  DATA  ACQUISITION  SYSTEM 

1.  Load  Real  Time  Executive  Basic  into  HP21MX  computer. 

2.  Load  KULITE  into  HP21MX.   Tune  pacer. 

3.  Put  disk  labeled  "Transonic  Compressor  -  Paige"  into  HP9867B  mass  memory. 

4.  At  start  of  a  run,  get  RESET1  from  unit  0  and  run  RESET1  to  initialize 
the  number  in  file  RECY//  to  unity. 

5.  Scratch  RESET1. 

6.  Get  TRAN4  from  unit  0. 

7.  Run  TRAN4.  The  HP9830A  display  will  remain  blank  while  the  HP9830A 
awaits  data  from  the  HP21MX. 

8.  Run  KULITE,  noting  that  the  two  mode  switches  on  the  Pacer  must  be  set 
according  to  the  A/D  converter  mode  to  be  used  (i.e.  0  or  4).   In  cali- 
bration, four  scans  of  the  twelve  channels  must  be  made.   The  first  scan 
must  be  with  the  pressure  at  S2  on  the  reference  side  of  the  Kulite 
transducers.   The  second  and  third  scans  are  made  with  other  steady 
calibration  pressures  on  the  reference  side.   The  fourth  scans  must  be 
made  but  any  signals  can  be  used  on  the  A/D  converter  provided  that  they 
do  not  cause  overloads.   In  mode  0  the  pacer  must  not  be  altered  during 
an  experiment. 

9.  On  completion  of  an  experiment  both  the  HP9830A  and  the  HP21MX  must  be 
reset  as  instructed  by  their  displays  prior  to  performing  another 
experiment. 

10.  On  completion  of  a  run  the  data  which  is  stored  temporarily  in  DATAY1 
file  must  be  transferred  to  permanent  files  CKRWm  duplicated  on  both 
unit  0  and  unit  1.   First  open  CRKWm  where  m  is  the  run  number.   The 
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length  of  CKRWm  must  be  set  at  k  records  where 
k  =  record  number  for  last  experiment  +  g 
and  g  =  1  if  last  experiment  is  a  calibration 

5   if  last  experiment  is  a  pacer  experiment 
13  if  last  experiment  is  a  free- run  experiment 
11.   To  abort  the  HP21MX  program,  enter  AB. 
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